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Professional Manpower for the War 


By HENRY A. BARTON 


Director, American Institute of Physics 


HE following outline is a study of the job a 
national manpower administration faces in 
a highly technical profession. It is idealistic—at 
least it probably goes beyond the hope of practi- 
cal realization. Nevertheless it was prepared as a 
contribution to the work of the branch of the new 
National Manpower Commission which is charged 
with administration of professions. 


I. ADMINISTRATION OF MANPOWER 


Nothing so important is so scarce as manpower 
of high professional caliber in certain fields whose 
utilization is essential to the. winning of the war. 
A national administration is imperatively needed 
to insure that available manpower in these fields 
is fully used and most effectively used. It is 
needed further to take all possible steps to 
augment the supply. Since positive measures in 
so complex a structure as a nation cannot be 100 
percent effective, the administration must also 
provide safeguards against known forces tending 
to waste the supply or use it at less than full 
efficiency. 


II. UNDERSTANDING THE PROBLEM 


Such a manpower administration can be suc- 
cessful only to the extent that it has an intimate, 
competent, and comprehensive understanding of 
the manpower situation in all of its aspects in 
those fields which are important. This implies 
more than statistics, however valid, complete and 
comparable. It is difficult to delimit what is 


implied. Included must be a real knowledge of the 
subject matter of the field, the pedagogical ex- 
perience in the field, the points of view and habits 
of thought of members of the profession, their 
organizations, the experietice of industrial man- 
agement in applying the field, the nature of 
services the field renders to offensive and de- 
fensive combat operations, etc. 


Ill. HOW TO OBTAIN THE UNDERSTANDING 
—EXAMPLE: PHYSICS 


Such an understanding can be pictured best by 
an example. It is proposed to outline below what 


‘is needed to obtain an understanding of a 


particularly compact, unified, and strategic field 
—namely, physics—with the thought that the 
new Manpower Commission should attempt to 
obtain this understanding, and, in so doing, 
establish a pattern of procedure, obtain experi- 
ence in fact-finding, and set up standards for the 
critical examination of information gathered—all 
as a valuable guide to the process of obtaining a 
similar understanding of other strategic fields. 


IV. QUANTITY, QUALITY, AND LOCATION 
OF SUPPLY OF PHYSICISTS 


This function has been developed by the 
National Roster and needs little discussion. A 
manpower administration may need more fre- 
quent reports as to location than have been 
available. One secretary, full time, could keep a 
list current if several other sources of information 
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(American Institute of Physics lists, government 
rolls, etc.) were coordinated with the Roster. 
This job can be done better and with less 
annoyance in the manner of a secretary keeping 
track of members than of a government agency 
using extensive questionnaire procedures. It is 
pertinent to enquire what records of physicists in 
uniform the Army and Navy maintain. 


V. PRESENT USE OF PHYSICISTS 


With little or no national coordination and 
insufficient understanding of the national supply 
and demand, the work of physicists is at present 
administered by various agencies. These include: 


A. The Army, Navy, and other armed services. 

B. Special government agencies coordinating research 
and development on war devices 
(O.S.R.D., N.A.C.A., N.R.C., etc.). 

C. Other government agencies supporting the war effort 
(National Bureau of Standards, Public Health Ser- 
vice, etc.) 


and methods 


1). Industries engaged in development and production 
of war devices and supplies. 

E. Industries making instruments and equipment neces- 

sary to research and testing work of A, B, C, and D. 

Universities, colleges, and high schools. 

G. Miscellaneous organizations including endowed re- 
search laboratories, national societies, independent 
consulting laboratories. 


— 


1. Individuals working as consultants, writers, etc. 


VI. EFFICIENCY OF USE 


The efficient use of a physicist requires the use 
by him of his intellect and energy in work which 
by its nature demands for its successful ac- 
complishment the knowledge, habits of thought, 
resourcefulness, and other qualities which are 
specifically developed by training and experience 
in physics. Such work is usually teaching, re- 
search, development, or operations but it may 
be administrative. Insofar as it is possible or 
desirable to distinguish the capacities of a 
physicist from those of an engineer, chemist, 
geophysicist, biophysicist, meteorologist, etc., a 
physicist has relatively more readily available for 
-his use a knowledge and understanding of other 
branches of physics than the one most used by 
him and an understanding of his own specialty 
which goes deep into the pertinent fundamental 
principles and facts and the mathematical re- 
lationships between them. 

A manpower administration would be obliged 


416 


to assure itself that the organizations of Section V 
were using physicists efficiently. This involves 
knowing, among other things: 


A. That every physicist in the organization is working 
as a physicist (at least any other use would require 
defense on the ground that it is either intrinsically 
more important work or is in a field in which the 
individual is at least equally competerit, and in which 
there is at least as serious a shortage of supply). 

B. That provision is made of adequate equipment. 

C. That where any portion of the work on which 
physicists are engaged can be done by others less 
trained in physics, such helpers are provided (labora- 
tory assistants, set wirers, mechanics, glass blowers, 
data takers, purchasing clerks, laboratory managers, 
etc.). 

D. That there is a wise administrative choice of what 
shall be taught, developed, and operated or in what 
directions research carried out—i.e., that physicists 
be set to the most useful tasks. Here again some coor- 
dination on a national scale is essential. Examples: 
(a) A physicist teaching a course.in a branch of 
physics not important to the war may be needed to 
operate a device to detect submarines. (b) A physicist 
who is “putting-together’’ radio-locaters may be 
more useful training hundreds of men to operate 
them. 

E. That physicists will be disturbed as little as possible 
(as long as their work is in a valuable direction) by 
transfers or threats of impending transfers. (Short 
period deferments are disturbing; uncoordinated 
competition for one’s services is disturbing.) 


A valuable contribution to physicists’ morale and 
efficiency will be made if a manpower adminis- 
tration can bring about a situation in which every 
physicist can feel he is playing his proper role in 
the war effort, and it is generally recognized that 
his work is patriotic and important. 


VII. ALLOCATION 


The most obvious and a very difficult problem 
of manpower administration is the distribution of 
emphasis between major functions, i.e., the de- 
termination of priorities and the allocation of 
quotas. In physics the problem is well illustrated. 
What percentage of the national resource of 
physicists should be allocated to Army operations, 
Navy operations, research and development, or 
methods and devices of war, industrial production 
of war equipment and ‘supplies, college, or high 
school level teaching? 

An intimate knowledge of the situation as-a 
whole is essential in finding the answer. To- 
morrow’s answer will be different from today’s. 
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For example, last year in physics it was desirable 
to continue education to the most advanced 
levels and to initiate researches for new war 
devices. We are entering a stage in which more 
emphasis should be placed on the operation of 
devices which are now nearing completion and 
hence on the training of large numbers of men in 
the first principles of physics required for such 
operations. 

To follow such comprehensive developments 
and trends with intelligent allocations of quotas, 
broadly informed policy committees in the stra- 
tegic fields will be necessary. Physics has a War 
Policy Committee in the American Institute of 
Physics available for consultation and this Com- 
mittee is already concerning itself with such 
problems. 


VIII. TRAINING 


The supply of physicists being insufficient, a 
manpower administration would coordinate facil- 
ities for training men. Experience with functions 
already discussed would indicate what aspects of 
physics should be taught (other than funda- 
mentals and prerequisite mathematics) and how 
far training should be carried. The number 
needed at various levels would have to be de- 
termined. The administration must see that 
teachers are used effectively in classes neither too 
small nor too large and that proper equipment 
for teaching and laboratory experience is avail- 
able (involving recommendations to the War 
Production Board regarding priorities). Above 
all the number of physicists needed for the 
teaching work must be determined and the con- 
ditions for their efficient use (see Section VI) 
satisfied (especially freedom from transfers— 
VI-E). 

Parenthetically there is a distinction between 
‘“‘a man who has had a basic course in physics” 
and a “‘physicist.’’ By the latter we mean a holder 
of the bachelor’s degree with major emphasis on 
physics plus one year or more experience or 
further education. The former may be trained for 
many operations, relieving some physicists. It is 
to be noted that short intensive educational 
procedures do not create a physicist except out of 
a man who is already ‘“‘almost”’ a physicist. 

The administration would require current 
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figures on enrolments of students and would 
probably need to assemble and distribute voca- 
tional information which schools, colleges, and 
universities could use to stimulate enrolments up 
to quotas desired. 

It is possible that an efficient national program 
would require the concentration of students, 
teachers, and equipment in a suitable number of 
selected institutions. Data enabling a wise selec- 
tion would have to be obtained. 


IX. INFORMATION AND GUIDANCE 


In addition to directions, the administration 
would have to issue information and guidance to 
organizations and individuals. Only from this 
source could the latter get information which is 
based on a national viewpoint. The morale of 
individuals can be improved by proper attention 
to their inquiries as to the best direction of their 
efforts. Information assembled by the adminis- 
tration for its other functions might suffice for 
this one. 


X. TECHNIQUE OF GATHERING 
INFORMATION 


Theoretically the manpower administration 
would need intimate information about the work 
of organizations, prying (e.g.) into the assign- 
ment of physicists to certain types of service in 
the army or into the respective merits of two 
projects of the National Defense Research Com- 
mittee. Practically this would amount to an 
absurd usurpation and centralization of the 
management function. It should be possible to 
establish, short of this, an understanding relation- 
ship with operating agencies in which any 
seeming interference with the latter would be 
more than counterbalanced by the benefits 
arising from the existence of a manpower ad- 
ministration. Such a degree of cooperation with 
these agencies can only be obtained if they are 
themselves really represented in making the 
judgments, decisions, allocations, and the like on 
the basis of the information assembled. The 
administration’s contact men should be compe- 
tent and acceptable to all agenciés. It is suggested 
that policy committees by fields should be called 
upon for the initial contacts as well as subsequent 
consultation. 
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 fgenrranpie is perhaps the most interesting 
and practical phenomenon in the theory of 
forced vibrations. The physicist and engineer are 
vitally concerned over the conditions under which 
resonance may occur, either for the purpose of 
eliminating it or encouraging it. For systems 
having one degree of freedom, three types of 
resonance are possible, namely: harmonic, mul- 
tiple harmonic, and sub-harmonic. These derive 
their names from the relationship between the 
frequencies of the applied force and the system; 
the frequency of the system being, respectively, 








Fic. 1. A pendulum with a horizontally moving (harmonic) 
point of support. 


equal to, an integral multiple of, a sub-multiple 
of, the frequency of the applied force. It is the 
purpose of this paper to discuss the mathematical 
and physical conditions under which each of these 
types of resonance may occur, and to illustrate 
each type with photographs of actual models. 


I. THEORY 


The equation of motion of a-vibrating system 
having one degree of freedom can be written in 
the form, 


md?*x /dt? = F,+ Fi+F., (1) 


in which m is the mass, x is the displacement, F, 
is the spring or elastic force, Fz is the damping 
force, and F, is the external or applied force. An 


418 





investigation of resonance involves a study of the 
response of the system described by the equation 


md?x /dt? = F,+ Fa, 


to the external force F.,. 


Harmonic Resonance! 
Consider Eq. (1) in the form, 


md?x /dt?+-Bdx/dt+kx = Fo sin wt. (2) 
This is a linear equation and its solution for large 
values of ¢ is easily found to be: 


x=A sin (wt—¢), 


in which the phase angle, ¢, equals tan! Bw/ 
(k—mw*) and the amplitude A_ equals 
Fy/((k—mw*)?+B*w? }'. If B?<2km, A is a maxi- 
mum when 


w?=wo?— B?/2m*, where wo?=k/m. 


(3) 

When this condition is satisfied the external 
force, Fo sin wt, causes our system to have its 
maximum amplitude, and we have the phenom- 
enon of harmonic resonance. If 8 is small we see 
that harmonic resonance occurs when the fre- 
quency of vibration of the external force is ap- 
proximately the frequency of the free vibration 
of the system. 


Multiple Harmonic Resonance’ 
Consider next Eq. (1) in the form, 


md*x /dt? + Bdx/dt+kx 


=)>A, sin nwt+ >-B, cos nwt. (4) 


This linear equation has, for large values of t, the 
solution : 


i sin (nwt — Yn) 
[ (k—mn*w?)?+ B?2n*w? |} 
B, cos (nwt — an) 


L(k—mn?w*)?+8?n*w? }} 








from which we see, provided 8 is small, that the 
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nth term in the sine and cosine series will have a 
very large amplitude if, 


Na" = wo” — B?/2m*. (5) 


Thus if our external force F, has in its expansion 
a term whose frequency, nw/2z, is related to the 
free vibration frequency by the above expression, 
resonance will occur. Therefore an ‘‘impure”’ 
applied force of frequency f can be in resonance 
with a system whose frequency is an integral 
multiple of f. This is the case of multiple reso- 


nance, or resonance at the higher harmonics. 


Sub-Harmonic Resonance 


The two previous cases were represented by 
linear* differential equations of motion in which 
the spring force was proportional to the displace- 
ment. This is no longer true for the case of sub- 
harmonic resonance. The special case of half- 
harmonic resonance is often observed in a system 
whose spring force is an explicit function of time; 
sub-harmonic resonance involving the sub-mul- 
tiples 3, }, 3; --- is observed in systems whose 
equations of motion are non-linear. 


Consider Eq. (1) in the linear form,? 
md*x /dt? + Bdx/dt+(k—2A sin 2wt)x=0, (6) 


in which m, 8, k, and Ao are constants. Knowing 
from experiments (von Melde’s, etc.) that a 
steady state solution exists, we can assume, if 8 


0 o ? f°? p Pp: iz 
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Fic. 2. Harmonic resonance. 
* An equation in the form d?x/d?@+ Pdx/dt+Qx=R is 


said to be linear provided P, Q, and R are functions of ¢ 
or constants. 
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Fic. 3. A pendulum maintained in motion by impulses 
from a revolving crank arm. 











and A, are very small, that the solution, to a first 
approximation, is in the form: 


x=A sin w+B cos ot. 
Putting this value of x in Eq. (6) it becomes, 


dD, sin wt—A oB sin Sut 
+D, cos wt+A oA cos 3wt=0, 


in which D,=kA—mw*A—BwB-—AjoB, and 
D.=kB—mw*B+8wA —AvA; and must be true 
for all values of t. Letting wt have the values 0 and 
1/2, we get the following equations: 


BwA +(k—mw*)B=0, (7) 
(k—mw*)A — BwB=0. 


The existence of our assumed solution depends on 
the existence of A and B, and therefore on the 
vanishing of the determinant of Eqs. (7). We find 
the condition for the existence of our solution to 
be: 


w? = B? 2m?+[(82/2m?)?+k?2/m? }. (8) 


For small values of 8 this reduces to the equation 
w*2w,”. Thus the applied force, (2A sin 2wt)x, 
can keep the system in a steady state of motion 
if the frequency (w/z) of the applied force is 
twice the frequency (wo/27) of the system. We 
have, therefore, an example of half-harmonic 
resonance in a linear system. 
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Consider next Eq. (1) in the form,’ 


md*x/dt?+-kx = F,(t) + f(x, dx/dt), (9) 
in which |e| <1, —kx+ef(x, dx/dt) = F,+ Fa, 


and k and e are constants. This equation is in 
general not linear, but it can be solved approxi- 
mately by the perturbation method. To this end, 
we expand F,(t) in a Fourier series, 


F(t) = >A, sin not + > B, cos not ; 
and assume that the solution of Eq. (9) is, 
x= xO 4 ex + Cx + Bx + - cn oe 


Putting this value of x in Eq. (9) and equating 
like powers of ¢, we have: 


md?*x /dt?+kx = >A, sin nwt+ > B, cos nut 
md?*x\ /dt?+kx =f (x, dx /dt) 


The first of these equations has as its solution, 
x — C sin (wot —a) 
A, sin nwt 


+— 


B,, cos nwt 
m(wo" — n*w") 


’ 


mM (wo? — n*w") 


in which wo? =k/m, and C® and a are constants 
yet to be determined. We now have x and hence 
dx /dt as a function of ¢, and thus our second 
equation can be written, 


md*x /dt?+kx™ = Ya, sin nwt+ > b, cos nwt, 
which has as its solution, 
x%=C™ sin (w°t—a) 

a, sin Nwt 


+5—————+5 


Mm (wo? — n*w") 


b,, COS Nwt 


M(wo" — n*w") 


The constants C™ and a‘” are chosen so that at 
any time ¢=0, the system is unperturbed; i.e., 
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Fic. 4. The discontinuous force F supplied by the 
motor in Fig. 3, as a function of time ¢. 
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Fic. 5. Multiple harmonic resonance. 











x) =dx /dt=0 when t=0. Thus to the first 
approximation, the solution of Eq. (9) is, 


x=x™ (CO, a ) + ex, 


in which C® and a are determined by the 
initial conditions on x and dx/dt. By continuing 
the above perturbation process, Eq. (9) can be 
solved to any degree of accuracy, the method 
being applicable provided ¢f;(x,dx/dt) can be 
considered as a term which ‘“‘perturbs’”’ the 
equation, 


md*x /dt?+kx = F,(t). 


To insure resonance in the first approximation, 
we must require x to have a frequency w/2r. 
Since x“ can be written in the form, 


t 
xY=(1 ma) f fila, dx™ /dt) sin wo(t—1)dr; 
0 


the conditions on C™ and a so that x“ does 
havea frequency w/27are given by the equations: 
~ 27 /wo 


filx(7), 


x) (294 wo) — 


0 


d{ x (r) ]/dt} sin wordr=0, 
(10) 


d[ x) (2x ‘wo) ] a= { fifx(r), 


0 


d[x(r) ]/dt} cos wordr=0. 


JOURNAL OF APPLIED PHYSICS 











ig 


0) 


CS 








The solution of Eq. (9) toa first approximation is 
now, 


x=x + ex (CM, gi), 


in which x“ has a frequency wo/2m7 and C™ and 
a” are to be determined by the initial conditions. 

If F.(t) is a force of small amplitude but fre- 
quency Nw»/2m, then approximately, x has a 
frequency wo/2m; but we made x“ (by Eqs. (10)) 
have a frequency wo/2z, and therefore our solu- 
tion x has a frequency wo/27. Thus the terms 
introduced by the damping and non-linear part 
of our equation, i.e., by «#0, provided real 





et 
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Fic. 6. A pendulum with a vertically moving 
(harmonic) point of support. 


values of C® and a® can be found to satisfy 
Eqs. (10), do not cause the natural vibration, of 
frequency w»/27 to vanish and be replaced by the 
forced frequency mwo/2m as would be the case if 
the system were linear. Therefore a non-linear 
system of frequency f is kept in motion by an 
external force of frequency nf. 


Il. EXPERIMENTAL EXAMPLES 


As an experimental example of harmonic reso- 
nance consider the system shown in Fig. 1, whose 
equation of motion is: 


d*6/dt?+(1/L)d*s/dt? cos 6+(mg/2L) sin 6=0. 


The massive vertical cantilever supplies approxi- 
mately simple harmonic motion s=A sin wt to 
the supporting end of the pendulum, so that the 
equation of motion can be written, 


d*6/dt?+K, sin @= Kz cos @ sin wt, 


in which K,=mg/2L and K2=Aw?/L. For small 
values of @ this becomes, 


d*?6/dt?+K,0=Kz sin ot, 
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and is the same as Eq. (2) with @=x, and B=0. 
Thus resonance will occur if the frequency of the 
pendulum is the same as the frequency of the 
cantilever. These results are shown in Fig. 2, 
which is made up of rapid sequence photographs‘ 
of system 1 in motion. 

The pendulum can be seen oscillating with 
considerable amplitude due to the small hori- 
zontal motion of its pivoted end. From pictures 
1, 13, 24, it is quite obvious that the pendulum 
and the cantilever have the same frequency. 

As an example of multiple resonance, consider 
the system (shown in Fig. 3) in which a pendulum 
is supplied with impulses having a frequency one- 
half the frequency of the pendulum. The equation 
of motion for small values of @ is, 


d*6/dt? +-w?0 = >°b, sin (nwt/T) O<t<T, 


in which the right-hand term represents the 
periodically applied force f(t) which is shown 
graphically in Fig. 4. After m (1, 2, 3, ---) cycles 
of this applied force the equation of motion 
becomes, 


d?6/dt?+wo?6 = >-b, sin (n/4m)wot ; 
in which the Fourier coefficients are given by the 


equation, 


4m /wo 
b,= (wo/2me) f f(t) sin woldt. 
0 


From our theoretical discussion of multiple reso- 
nance we know that it will occur if there is a term 
in f(t) which has the same frequency (wo/27) as 
has the pendulum. This is so when n=4m; and 
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Fic. 7. Half harmonic resonance. 
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show 


to guarantee resonance we must that 
bam #0. Putting the values of f(t) in the equation 
for b, we find: 


(4rj/wo)+e 


b,= Fo dD 


j=1 ixj/wo 


sin wol dt 
= (Fo /2n)(1 — COs €w). 


Thus 54,40, and therefore we have an example 
of multiple harmonic resonance, the results of 
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Fic. 8. A cantilever spring of gradually increasing stiff- 
ness maintained in motion by an unbalanced motor 
mounted on its free end. 


which are shown in Fig. 5. From pictures 1, 13, 25, 
it is quite obvious that the applied force has a 
frequency one half that of the pendulum. 

A system exhibiting half harmonic resonance 
is shown in Fig. 6. Its equation of motion, for 
small values of 8, and a massive pendulum, is 


d*6/dt?+- (wo? —K sin wt)@=0, 


in which wo?= g/L and K = —Aw?/L, L being the 
length of the pendulum, A» the maximum ampli- 
tude of the applied force, and w its angular 
velocity. This is the same as Eq. (6) with x=@ 
and 8=0, and we expect resonance if w is twice 
wo. These results are shown in Fig. 7. Pictures 
1, 11, 21, show that the applied force has twice 
the frequency of the pendulum. 

An example of sub-harmonic resonance,’ shown 
in Fig. 8, has as its equation of motion, 


(M+m)d*z/dt?+[a—(M+m)g/R ]z 


= —mrw* cos wl — Kz’, 


in which z=R8, w is the angular velocity with 
which m is rotated about M, and —az— Kz’ is the 
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non-linear spring force resolved into a straight 
line and a cubical parabola. The actual spring 
force is shown by curve I, Fig. 9; the linear and 
non-linear components are shown by curves II 
and III. For this model, a= 33,300 dynes/cm and 
K=1700 dynes/cm*. The above equation of 
motion can be written in the form, 


d*z/dt?+ (wo)? = B cos wt — ez’, 


in which wo’ =[a/(M+m)-—g/R}, B= —mrw*/ 
(M+m) and e= K/a. This is the same as Eq. (9) 
and has as its zeroth order solution, 


2 =C sin (wo/t—a) 


+[B/[(wo’)? —w?]] cos wt, 


in which C® and a must be determined so that 
z has a frequency wo /27. If we consider the case 
of w=3wo’, the equations determining C® and 
a® are: 


2r/wo’ 
f [co sin (wot — a) 
+B’ cos 3wo't |® sin wo'tdt=0, 


27 /wo’ 
f [C® sin (wo’t—a) 


0 
+B’ cos 3wo't |* cos wo’ tdt=0, 
in which B’=B/[(wo’)?—w? ]. From these inte- 


grals we determine that real values of C™ and 
a exist. Thus we have an experimental example 


F 











FyG. 9. Curve I. The non-linear spring force F, of the 
cantilever of Fig. 5, as a function of displacement z. This 


is resolved into a linear curve (II) plus a non-linear curve 
(III). 
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Fic. 10a. Sub-harmonic resonance, one to three. 


of one-third harmonic resonance. The experi- 
mental results are shown in Fig. 10a. Pictures 
1, 4, 8, 12 show that m makes three revolutions 
while M completes one cycle. 
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Fic. 10b. Sub-harmonic resonance, one to four. 


The experimental results for (}) sub-harmonic 
resonance are shown in Fig. 10b. Pictures 1, 3, 6, 
8, 10 show that the applied force has four times 
the frequency of the cantilever. 
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Résumés of Recent Research 








New Studies on the 
Clusters of Nebulae 


Investigations on the 
large scale distribution 
of matter in the universe 
were concerned until recently with the geometrical 
or geographical aspects of this distribution, 
which, chiefly through the work of Hubble!’ and 
others, can now be considered as established in 
general to distances of some hundred millions of 
light years. The problem of current interest is the 
analysis of the physical principles governing the 
distribution of nebulae. For example, the question 
may be raised if Newton’s law of gravitational 
attraction still holds good when the interacting 
masses are separated by distances of millions of 
light years rather than by distances not ex- 
ceeding a small fraction of one light year such as 
are encountered among the members of the solar 
system and the components of those double stars 
for which the validity of the law has been checked. 

Observations on the clusters of nebulae prom- 
ise to furnish proof for a true universality of New- 
ton’s law to an excellent degree of approxima- 
tion. These observations refer to clusters which 
have reached statistically stationary states; i.e., 
clusters satisfying the following conditions: 

(a) The clusters must possess spherical sym- 
metry subject only to fluctuations derivable 
from the limited number of nebulae involved. 

(b) The velocity dispersion throughout the 
cluster must be independent of the distance from 
the center out to distances where the mean free 
paths of the nebulae exceed a certain critical 
length (Smoluchowski limit). 

(c) The radial distribution curves of the num- 
bers of nebulae per unit volume for different 
clusters must be reducible to the same standard 
Emden curve which represents the distribution of 
matter in a bounded isothermal gravitational gas 
sphere. The reduction factors involved must be 
derivable from the velocity dispersion and the 
central density of the clusters. 

(d) A segregation of nebulae of different mass 
must occur. 

Available data on the clusters in Coma, Hydra, 
Perseus, Cancer, Fornax, and Pegasus indicate 


1 E. Hubble, The Realm of the Nebulae (Yale University 
Press, 1936). 
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that these clusters satisfy the mentioned require- 
ments and that they are stationary assemblies of 
nebulae whose interactions are governed by the 
inverse square law of gravitation. 

Analysis of the observations on clusters of 
nebulae also indicates that the time scale neces- 
sary for a rational understanding, on the basis of 
fundamental principles of physics, of the ob- 
served large scale distribution of matter is much 
longer than the hypothesis of the expanding uni- 
verse would have it.’ 

2 F. Zwicky, Phvs. Rev. 61, 489 (1942); Astrophys. J. 86, 
217 (1937), and May issue, 1942. Proc. Nat. Acad. Sci. 25, 
604 (1939); 26, 332 (1940); 27, 264 (1941); 27, 366 (1941); 


see also April issue, 1942. Publ. Astr. Soc. Pacific 50, 218 
(1938). 


New Micro- 
photometer 


The introduction into 
industrial service of pro- 
cedures developed in the 
laboratory usually entails extensive modification 





Fic. 1. \ production-type microphotometer. 
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Fic. 2. The optical system. 
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of apparatus design to provide the ruggedness 
and stability required for shop use. In these re- 
spects, spectrochemistry has proved no excep- 
tion; current improvements in instrumentation 
for specific applications are numerous. The new 
microphotometer described! by Vincent and 
Sawyer is adapted especially to installations 
where durability and continuous reliability are of 
more importance than initial cost. Specifically, it 
is an attempt to provide an instrument adequate 
for routine ferrous laboratories which may be 
called upon to make a million readings per year at 
the highest attainable precision. 

The instrument comprises a two hundred 
pound cast base carrying rails which support a 
stage adapted for motion along the spectrum. 
The stage, in turn, carries on a vee way and steel 
ball a cross-slide on which is mounted a spring- 
loaded clamp for the spectrum plate. The line- 
arity of the longitudinal motion is provided by a 
pyramidal rail and two grooved rollers in front; 
the rear support is a single cylindrical roller and 
flat-topped rail. Rapid positioning of a spectral 
line for measurement is made possible through 


1H. B. Vincent and R. A. Sawyer, “‘A new microphotom- 
eter,”’ J. Opt. Soc. Am. 31, 639-643 (1941). 


the use of two frictionally loaded work wheels 
attached to the stage and engaging a worm 
rotably mounted between two supports attached 
to the base. The worm may be rotated by either 
of two coaxial hand wheels, shown in Fig. 1, con- 
nected to it by gear trains of different ratios, 
providing, in all, three speeds of longitudinal 
motion. Gears are spiral, quarter-inch face, 
hardened and polished; shafts and rollers are 
carried on ball bearings. 

Excessive error from emulsion grain is avoided 
by photometering a line length of five millimeters 
at as great a width as is consistent with resolution 
requirements. Actual operating widths vary from 
35 to 50 microns. The optical system, deriving its 
illumination from a thirty-watt lamp, is shown in 
Fig. 2, which details also the provision made for 
the projection of a reference spectrum carried on 
the stage. The indicating instrument is a galva- 
nometer of short period located at two meters 
from its scale to reduce obliquity errors. Close 
grouping of projected images and scale are de- 
signed to minimize fatigue. 

A test on the reproducibility of reading of a 
line pair indicates a relative error of less than 
0.05 percent in the average ferrous analysis. 





Here and There 


Honor Roll 








Five scientists were lost in a collision of two blimps in an 
accident off Manasquan, New Jersey, Monday night, 
June 8. The men were: L. S. Mover, Assistant Professor of 
Botany, University of Minnesota; C. R. Hoover, Pro- 
fessor of Chemistry, Wesleyan University; F. C. GiLBert, 
New London, Connecticut; A. B. Wyse, San Diego, 
California; I. H. TrLLes, San Diego, California. In com- 
menting upon the tragedy, the Office of Scientific Research 
and Development has issued the following statement: 
“These men lost their lives in an accident encountered in 
an experimental flight in connection with important scien- 
tific war research and gave their lives for their country as 
truly as any soldier or sailor killed in battle.” 


* 
Honors and Awards 


Among those recently elected to membership at the 
annual raeeting of the National Academy of Sciences were 
the following: LYMAN J. BricGs, Director, National Bureau 
of Standards; ALBERT EINSTEIN, Professor of Mathematics, 
Institute for Advanced Study, Princeton; PAuL Dyer 
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Merica, Vice President, International Nickel Company; 
JoHN TorRENCE TATE, Dean of the College of Science, 
Literature, and Arts, University of Minnesota. 


The American Academy of Arts and Sciences has elected 
as fellows in the mathematical and physical sciences: 
WILMER LANIER BaARROw, the Massachusetts Institute of 
Technology; FRANCIS Bircu, Harvard University; SAMUEL 
CoRNETTE COLLins, the Massachusetts Institute of Tech- 
nology; Otto StRuVE, Director, Yerkes Observatory. 


At the annual meeting of the American Philosophical 
Society members elected in the mathematical and physical 
sciences were: OLIVER ELLSWORTH BUCKLEY, President, 
Bell Telephone Laboratories, Inc.; LEE ALVIN DUBRIDGE, 
Professor of Physics and Dean of the Faculty of Arts and 
Science, University of Rochester; DUNCAN ARTHUR 
MacInnes, Rockefeller Institute for Medical Research; 
ROBERT RAYNOLDS McMatn, Director, the McMath- 
Hulbert Observatory, University of Michigan; FRANCIS 
Dominic MURNAGHAN, Professor of Applied Mathe- 
matics, Johns Hopkins University; HAROLD MALCOLM 
WESTERGAARD, Gordon McKay Professor of Civil Engi- 
neering and Dean of the Graduate School of Engineering, 
Harvard University; ROBERT RUNNELS WILLIAMS, Chem- 
ical Director, Bell Telephone Laboratories, Inc. 
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Honorary members elected to the Academy of Sciences of 
the U.S.S.R. on May 8 included the following: Ernest O. 
LAWRENCE, Director of the Radiation Laboratory, Uni- 
versity of California and Gi_pert N. Lewis, Professor of 


Chemistry and Dean of the College of Chemistry, Uni- 
versity of California. 


The Optical Society of America has announced the 
award of the Adolph Lomb Medal for 1942 to JAMeEs G. 
BAKER of the Harvard College Observatory. This medal, 
which is awarded biennially for noteworthy contributions 
to optics, will be presented at the time of the annual meeting 
of the Society next October. The recipient has made im- 
portant contributions to the improvement of instruments 
by which astronomical observations are made. He has 
originated new methods of procedure in the fields of 
geometrical optics and lens design. 


Three young scientists have been awarded fellowships to 
carry on their investigations at the Westinghouse Research 
Laboratories, according to Dr. E. U. Condon, Associate 
Director of the Laboratories. 

The appointees are Dr. Russe_y E. Fox, who will carry 
on research on the dynamics of high speed machinery which 
he began at the University of Virginia under Professor J. W. 
Beams; DR GERSON S. SCHAFFEL, Carnegie Institute of 
Technology, who will engage in fundamental research on 
the chemistry of polymers; and Dr. Rospert S. WEIsz, 
Cornell University, who will work in the field of ceramics, 
attempting to improve the properties of electric porcelain. 

This is the fifth group to be selected under the Westing- 
house Research Fellowship plan which was inaugurated in 
1938. Appointments under the Fellowship plan are for one 
vear with reappointment for a second year if mutually 
desirable. Although in the past these men have devoted 
their entire time to investigations in the realm of pure 
science, this year they have been asked to indicate the 
fields in which they feel qualified to contribute to war re- 
search if their assistance should be required. Three of last 
vear’s appointees are now engaged in war research. 

Of last year’s fellowship winners, Dr. THomAs W. DAKIN 
Harvard University, has been reappointed for a second 
vear’s study. Dr. James S. KoenLer of the University of 
Pennsylvania has been granted a year’s leave of absence to 
act as instructor of physics at Carnegie Institute of 
Technology. The three other 1941 appointees have joined 
the Westinghouse Research staff. They are Dr. DANIEL 
ALPERT, Stanford University; Dr. JoHN W. COLTMAN, 
University of Illinois; and Dr. THEODORE HOLSTEIN, New 
York University. 


* 


Conference on Problems of Modern Physics 


To celebrate the public opening of the Charles Benedict 
Stuart Laboratory of Applied Physics, a Conference on 
Problems of Modern Physics was held at Purdue Uni- 
versity on June 19th and 20th. The program was as follows: 

JUNE 19 
Physics in Industry, Dr. E. U. CONDON, Associale 
Director, Westinghouse Research Laboratories. 


Opening address 
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JUNE 20 

Theories of Nuclear Forces, Dr. J. ScHwinGER, Purdue University. 

Problems of Modern Field Theories, Dr. W. PAautt, The Institute for 
Advanced Study. 

Design and Construction of the Twenty-Million-Volt Induction Ac- 
celerator, Dr. D. W. Kerst, University of Illinois. 

The Physics of Statistical Noise, Dk. W. W. HANSEN, Stanford Uni- 
versity and the Sperry Gyroscope Company. 

The Electron Microscope and Its Applications, Dr. J. A. 
Bell Telephone Laboratories. 


BECKER, 
* 


Summer Courses in Physics 


To help meet the present urgent need for people with 
advanced training in physics, the Graduate School of New 
York University, for the first time in a number of years, is 
offering courses in physics and mathematics during the 
summer. 

These courses will be given in two short terms from 
June 29 to August 7 and from August 10 to September 18. 
They are planned as part of a one-year curriculum designed 
to cover the basic subjects of physics, including dynamics, 
electricity and magnetism, thermodynamics, physical 
optics, atomic structure, advanced laboratory, and methods 
of mathematical physics. The course in thermodynamics 
will be given in the first summer term and that in advanced 
laboratory in the second. The course in methods of mathe- 


matical physics will continue through both terms. 
* 


From College Machine Shop to War Factory 


Purdue University’s student machine shops at Lafayette, 
Indiana, customarily used only for routine practice work, 
recently have been converted into a war production factory 
under the terms of a Westinghouse subcontract. 

More than 275 engineering students enrolled in the shop 
course are working part-time on the subcontract. Together 
their production is the equivalent of a 75-man machine 
shop working full time on war materials. Special stock 
material required for some parts being made by the 
students has been supplied by Westinghouse, and the uni- 
versity has received a priority rating from the War Pro- 
duction Board for the purchase of other materials. 


* 


Necrology 


Dr. MILLARD F. MANNING, Assistant Professor of 
Physics at the University of Pittsburgh, died June 1st of a 
ruptured appendix at the age of 36. For the past few years 
Dr. Manning had been interested in theoretical calculations 
concerning the properties of metals and other solids. Re- 
cently he and Dr. M. E. Bell completed an extended study 
of the properties of semi-conductors which was published in 
the Reviews of Modern Physics. Dr. Manning was a member 
of the American Physical Society. 
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Surface Replicas for Use in the Electron Microscope 


VINCENT J. SCHAEFER AND DAvip HARKER 
Research Laboratory, General Electric Company, Schenectady, New York 


(Received April 9, 1942) 


A simple method is described for makjng replicas of 
surfaces for study in the electron microscope. The surface 
to be studied is cleaned and then coated with a very thin 
slightly wedge-shaped film of polyvinyl formal having a 
range in thickness from about 500 to 750A in 2.5 cm. This 
film is formed by dipping the prepared surface into a 0.5- 
percent solution of the resin dissolved in dioxane or ethy- 
lene dichloride and permitting the solvent to evaporate 
while the surface is held in a vertical position. The replica 
film is removed by placing the coated surface under water 
and peeling the film from it. A description of this tech- 
nique is given in detail. A graph is presented which shows 
the thickness of resin films obtainable from various con- 


N a recent paper! Zworykin and Ramberg de- 

scribed a method for studying the surfaces of 
solid materials, such as the etched surfaces of 
metals. The authors pointed out that the methods 
described by H. Mahl** in a series of papers in 
1940 had the disadvantage that the specimen 
was destroyed in obtaining the replica. Zworykin 
and Ramberg described several techniques which 
could be used without spoiling the specimen and 
illustrated by micrographs the results obtained. 
Of the three methods described in some detail in 
their paper they indicated a decided preference to 
the silver-collodion process. This method involves 
the evaporation in a vacuum chamber of a metal 
film onto the surface of the specimen to be 
studied, which is then, if necessary, reinforced by 
electroplating, to such a thickness that the metal 
film can be stripped mechanically from the sur- 
face. A thin collodion film is then flowed over the 
metal replica, after which the latter is removed 
from the plastic by dissolving it in nitric acid. 
The collodion film is then mounted for study in 
the microscope. 

Several other methods were tried by them, in- 
cluding one in which a resin material termed 
“vinylite lacquer’’ was applied directly to the 


1V. K. Zworykin and E. G. Ramberg, J. App. Phys. 12, 
692 (1941). 

2H. Mahl, Zeits. f. tech. Physik 21, 17 (1940). 

3H. Mahl, Metallwirtshaft 19, 1082 (1940). 
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centrations of resin in dioxane, and a table gives the 
interference effects which may be used for determining 
approximate thicknesses of thin resin films. Several easy 
procedures are described for checking the thickness of 
such films. A technique is described for preparing strong 
films having a thickness of 100A and an area of 3 cm? 
when mounted in air. Such films are very useful when 
working with the electron microscope and several methods 
are described for mounting specimens on them, including 
a conditioning treatment—useful in mounting water 
suspensions—which makes one surface very hydrophilic. 
Five micrographs show the results obtainable with the 
techniques described. 


specimen and removed mechanically. This latter 
method was evaluated in their paper! by the 
following statement: ‘“‘While extremely simple it 
is relatively uncertain and yields far poorer 
resolution than the first method (silver-collodion) 
discussed.”’ 

This paper will attempt to point out possible 
reasons why this direct method can be “‘relatively 
uncertain” and possible ways in which this 
uncertainty can be nearly, if not completely, 
eliminated. It is obvious that maximum possible 
contrast can be obtained from a replica made 
directly from the original sample if a good replica 
is produced. Such a replica must be free of 
imperfections due to surface tension effects, thin 
enough to avoid electron absorption, thick 
enough to be removed without undue trouble, 
and strong enough to be removed from an etched 
surface without showing strains. 

We have found a material which seems to fit 
these requirements quite well, and Figs. 1-3 
illustrate the type of results which can be ob- 
tained. The method of forming such replicas is 
simple and the equipment needed can be as- 
sembled in a few hours. The resin used is polyvinyl 
formal known by the trade name of Formvar 
15-954 and is dissolved in dioxane or ethylene 





4 Shawinigan Products Corporation, Empire State Build- 
ing, New York City. 
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Fic. 1. Replica of etched surface of 0.4 percent C 
steel oil quenched. 


dichloride. This material has been used previ- 
ously® ® where very thin strong films were needed, 
and occupies an eminent place in certain indus- 
trial applications because of its outstanding re- 
sistance to abrasion and other mechanical stresses 
and strains. 

In preparing replicas of solid surfaces for obser- 
vation in the electron microscope certain limi- 
tations are imposed by the physical properties of 
the electron beam impinging on the specimen as 
well as other related factors. We have found that 
a film having a thickness of 1500A begins to 
deteriorate within a few minutes when used in 
an RCA type B electron microscope at 50-60 kv. 
With a film of this thickness, observation shows 
that the degree of contrast is quite poor. A study 
was made of the effect of thickness on contrast 
and the results obtained indicate that a thickness 
between 500 and 750A produced the best image 
formation with the power supply indicated above. 
Films whose average thickness was less than 300A 
showed contrast as poor as those whose thickness 
exceeded 1000A. To understand this anomaly it 


5 L. H. Germer, Phys. Rev. 56, 58 (1939). 
®V. J. Schaefer, J. Phys. Chem. 45, 681 (1941). 
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would be best to consider the factors involved in 
preparing good replicas. 

A surface carefully polished with MgO is rela- 
tively free of scratches. Those which remain are 
wide and shallow; the difference in height be- 
tween the deepest groove and the highest eleva- 
tion is probably not greater than 200A. From 
such a surface it is easy to obtain a stripping 
having a thickness of 250A. Under the electron 
beam little, if any, surface detail can be observed. 
With a light etch (2—5 sec. with a saturated solu- 
tion of picric acid in alcohol) on an oil-quenched 
steel sample containing 1.1 percent C, a stripping 
having a thickness of 300-350A can be obtained 
without difficulty. In the electron microscope 
such a stripping shows a certain amount of 
“background” indicative of the crystalline struc- 
ture of the metal, but the contrast is insufficient 
to be of any use to the observer. When the 
etching process is carried to the point where it is 
possible to get a fair photomicrograph at 1000X 
with the light microscope (15—25 sec. etch) it has 
been found that the replicas which can be ob- 
tained have a minimum thickness of about 550A. 
When an effort is made to secure replicas thinner 
than this latter dimension, the results are not 
successful. An increase of film thickness up to 
about 750A was found to be permissible without 
serious loss in contrast. These thicker films are 
also easier to remove from the metal surface, the 
main disadvantage being the loss in illumination 
on the fluorescent screen. It was found to be 
to work with films thicker than 
1000A because of this latter difficulty as well as 
the deterioration which occurs due to electron ab- 
sorption. Our results indicate that nothing is 
gained by working with replicas having a thick- 
ness greater than that necessary for obtaining a 
coherent film. The advantage in additional 
strength which is gained by having the replica 
thicker by 50-100A than the distance between 
the high and low points of the etched surface is 
probably a real one, but anything greater than 
this can only result in a loss of contrast. 

In order to have sufficient contrast to form a 
good image of an etched surface it is very im- 
portant that the upper surface of the resin film is 
flat while that in contact with the metal surface 
follows it in intimate detail. That surface tension 
forces actually produce such conditions can be 
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shown by preparing a relatively thick replica of a 
heavily etched surface of pure copper. A 3-percent 
solution of polyvinyl formal was used, the replica 
made in the usual manner, removed from the 
copper block, mounted on a frame and viewed in 
the light microscope at about 200X magnification 
with vertical illumination which produced inter- 
ference colors in the replica film. This film had an 
average thickness of about 3000A and under the 
microscope appeared identical to the original 
etched copper surface except that each grain 
exhibited certain interference colors. The im- 
portant fact noted was that the interference 
colors were quite uniform for a given crystal 
grain right up to the grain boundary. The ad- 
joining grain showed another color, also quite 
uniform over its entire area and was itself 
adjacent to other colored patches representing 
still other grains. A comparison of the observed 
interference colors with a thickness gauge made 
it an easy matter to ascertain the relative posi- 
tion in elevation of each grain but more im- 
portant for present considerations demonstrated 
that the upper surface of the replica was quite 
level and did not follow the ‘hills and valleys” 
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Fic. 2. Replica of etched surface of 1.1 percent C 
steel furnace cooled. 
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Fic. 3. Replica of etched surface of 0.4 percent C 
steel furnace cooled. 
of the metal surface. If this latter condition were 
obtained, the interference color for a given grain 
area would be considerably different in the 
center from that at the grain boundary. A change 
in thickness of less than 25A would be readily 
detectable under the conditions described. 

It is desirable from time to time to have a 
reasonably accurate and simple method for de- 
termining the thickness of the replicas which are 
being prepared from a stock solution so that its 
strength can be adjusted as the solvent disap- 
pears due to evaporation. The film thickness 
which our experience has shown produces replicas 
having the best contrast is made from a solution 
containing 500 mg of polyvinyl formal dissolved 
in 100 cc of dioxane. As shown in Fig. 4, the mid- 
thickness of a wedge-shaped film 2.5 cm wide 
made from a solution of this concentration, in a 
manner to be described later in- this paper, is of 
the order of 650A. The strength of the stock 
solution used may be checked occasionally by 
dipping a cleaned (with wet MgO on a cotton 
swab) and dried (with filtered air) 1X3” 
X 0.060” chromium plated slide’ into the solution 


7 Appollo Metal Works, 66th Place and South Oak Park 
Avenue, Chicago, Illinois. 
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and letting the solvent evaporate with the plane 
of the slide vertical and the long axis horizontal. 
If, when viewed in white light at nearly grazing 
incidence (i=80°) the thicker edge of the film 
coating the chromium plated slide appears as a 
yellow interference color, the solution has the 
proper concentration for producing good replicas. 
Table I represents the relation of observed color 
to actual thickness of a resin film flowed on a 
polished chromium surface. These colors are best 
observed by polarized light, although this refine- 
ment is not essential. From this table and Fig. 4 
it is apparent that the solution should not be 
allowed to become more concentrated than 0.75 
percent which would produce a film having a color 
range in a distance of 2.5 cm of yellow to purple 
when viewed at nearly grazing incidence i= 80°. 
The treatment of the surface of a solid prelimi- 
nary to forming its replica governs to a con- 
siderable degree the ease with which the replica 
may be removed. The replica of the surface of 
polished glass is probably about the easiest to 
obtain if the glass is reasonably clean. From such 
a surface it is a simple process to obtain a 
stripping having a thickness of 100A. A single 
layer of molecules of a hydrophobic substance, 
such as grease or oil, will completely prevent 
carrying out this normal stripping process and it 
is necessary to increase the film thickness at least 
tenfold in order to add sufficient tensile strength 
to the film to enable its mechanical removal. 
Metal surfaces, particularly those etched after 
polishing, present a much more difficult problem 
in preparing satisfactory replicas unless certain 
important factors are recognized. It was found, 
for example, that a delay of no more than five 
minutes between the drying of a steel specimen 
after etching and the formation of the resin film 
on its surface was sufficient to prevent the re- 
moval of the film from the metal. When this film 
was subsequently removed by a solvent and the 
surface of the metal cleaned by rubbing lightly 
with MgO on a cotton swab followed by washing 
_with distilled water and quick drying with a 
stream of filtered air, a similar film applied with- 
out delay was removed easily and showed ex- 
ceptionally fine detail. Great care must be exer- 
cised in the final cleaning of the metallographic 
specimen so that as the specimen is dried it is 
quite free of even a portion of a deposited 
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TABLE I. Actual thickness of thin films on a polished 
chromium plated surface having a refractive index 
range of n=Ca 1.49. 


Color at i =80° 


Actual thickness 





600A Pale yellow 
700A Yellow 
800A Yellow orange 
850A Orange 
900A Orange red 
950A Red 
1000A Purple 
1050A Deep blue 
1100A Blue 
1150A Light blue 
1200A Pale blue 
Faint blue 


1300A 


monolayer of hydrophobic types of molecules. It 
is very difficult to obtain this condition with the 
usual types of thermoplastic mounts often used 
for holding the specimen during the polishing 
process since the bond between the metal and the 
mounting medium often pulls away to leave a 
tiny crack which harbors surface active mole- 
cules. As the specimen is dried these molecules 
tend to spread over the metal to form metallic 
soaps and other types of adsorbed monolayers 
which later are responsible for lack of success in 
obtaining good replicas. These layers often will 
form in less than a second. A good test for the 
presence of a contaminating or hydrophobic layer 
on the metal surface is to place a drop of distilled 
water from a carefully cleaned pipette on the 
surface. If it flashes out and covers the surface 
without a contact angle, the surface treatment 
may be considered adequate. If, however, the drop 
remains as such and shows a definite contact angle 
with the metal, the cleaning method must be im- 
proved upon if good replicas are to be obtained. 
For preparing good replicas the following tech- 
nique has been found to be very successful : 
After the specimen has been given a good 
metallographic polish (we use MgO on polishing 
cloth known as kitten’s-ear for the final polishing 
of steel specimens), the surface is given an etch of 
such a depth that a fairly good photomicrograph 
can be obtained at 1000X in the light microscope. 
The etching solution is then removed and the 
surface quickly dried with a jet of clean filtered 
air. The specimen is then immediately lowered 
face down into a petri dish containing a small 
amount of a 0.5-percent solution of polyvinyl 
formal dissolved in dioxane. The specimen is 
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removed with a smooth motion and swung into 
a position such that the plane of the surface is 
vertical with respect to the surface of the resin 
solution. The solvent is allowed to evaporate to 
completion while held in the vertical position. 
This procedure forms a wedge-shaped film on the 
metal surface thickest at the lower edge and 
thinnest at the top. With a specimen having a 
diameter of 2.5 cm the range in thickness will be 
from 500 to 750A, although at the very bottom 
there will be a small area much thicker than this 
where the draining liquid accumulates. The 
specimen bearing the resin film is then placed 
under water in a shallow tray at such a depth that 
the metal surface is about 3 cm under the water 
surface. 

A good pair of sharp pointed tweezers with 
both “‘anvils’” in good register is then brought 
into contact with the edge of the metallographic 
specimen at the thick end of the resin wedge and 
the film is carefully “‘teased’”’ until it is seen to 
pile up in a small roll. After this has proceeded 
for a millimeter or so toward the center of the 
specimen, a purchase is obtained on the edge of 
the film and the tweezers are slowly drawn to- 
ward the thin end of the wedge. 


THICKNESS OF WEOGE-SHAPED RESIN FILMS 
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Fic. 4. Relation of concentration of resin in solvent to 
thickness of wedge-shaped film. 
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If the surface was free of hydrophobic material, 
the etch of the right depth, and the metal rela- 

tively free of porous inclusions such as oxides, a 

considerable piece of replica film will easily peel 

from the surface. Sometimes the complete film 

may be removed, although as a rule the best 

replicas tear before reaching the opposite side. 

While still under water, the peeled film is swung 

around so that the upper surface will reach the 

water surface first in an extended condition. Just 

before bringing the film to the surface, the latter 

should be swept clean of any contaminating sur- 

face film. As the film arrives at the surface it will 

shed water completely and should be entirely 

free of wrinkles or folds. While still on the water 

surface, the nickel mesh screens for mounting the 

replica in the microscope may be positioned in 

selected positions on the floating film, after which 

the film is removed from the water surface. In 

doing this, we prefer to place a bit of adhesive 

material, such as rubber cement, on the rim of a 

metal frame which is then brought down onto the 

upper surface of the replica. The frame bearing 

the film and screens is then removed from the 

cleaned water surface with a smooth tipping 

motion which finds the plane of the resin film at 

right angles with the water surface at the end of ‘ 
the movement. This procedure removes the film 

from the water without subjecting it to undue 

strain and in a dry condition. With certain 

metals, such as duralumin, where one constituent 

is removed preferentially by a special etching 

technique,® the residue will be removed with the 
replica film and consequently the under surface 
of the film will be wet when removed from the 
water bath. A convenient set-up for carrying out 
these steps is a 6” petri dish 2}’’ deep painted 
black on the outside and illuminated by the 
reflection of white light from a sheet of white 
paper or cloth. The edge of the dish should be 
made hydrophobic by rubbing with molten wax, 
and two narrow strips of brass 3” X}/’ X8” simi- 
larly treated will be useful for cleaning the water 
surface of surface contamination. 

In the course of making replicas of metal sur- 
faces with polyvinyl formal dissolved in ethylene 
dichloride, a peculiar spotted structure was often 
observed. At first it was thought that the ob- 





8F. Keller and R. A. Bossert, Metal Progress 41, 63 
(1942). 
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served structure was characteristic of the metal. 
This idea was discarded, however, when similar 
detail was found on various unrelated metals in- 
cluding pure nickel and copper. Other reasons 
were suggested to explain the occurrence of the 
spots which sometimes varied in size and density, 
but each was in turn successively discarded. 
Finally, a theory was proposed that the spots 
might be due to an oleophobic property of the 
surface which by a mechanism previously studied® 
would cause the liquid solution of resin to start 
peeling back from the metal at many nuclei. Be- 
fore the film could thin sufficiently to form actual 
holes, as occurred in the work cited, the solvent 
evaporated. 

In order to test out the validity of this pro- 
posed explanation of the spotted structure, a 
chromium plated slide well polished was selected 
and given. various surface treatments. When 
cleaned very carefully with MgO, taking every 
precaution to prevent the deposition of a mono- 
layer, a replica from the surface was nearly free 
of spots and those present were very small. A 
monolayer consisting of melted ferric stearate 








Fic. 5. Particle of stearic acid smoke mounted between 
two thin films of polyvinyl formal. 
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well rubbed was then tried and this time the 
replica showed the type of spotted structure 
which was typical of many previous strippings. A 
very oleophobic® surface was then prepared of 
two monolayers of acid barium copper stearate. 
The replica from this surface showed unusually 
large spots again spaced quite at random. 

Although the presence of spots proved to be 
undesirable in the present instance, it should be 
mentioned at this point that it will be worthwhile 
to study this effect more carefully since it is quite 
possible that the presence of the holes, especially 
when they show a definite alignment, might be 
showing the location of oxide or other types of 
inclusions so small as to be beyond the resolving 
power of the present electron microscope. 

By using dioxane instead of ethylene dichloride 
as a solvent for polyvinyl formal resin, the 
trouble due to spotting was entirely eliminated 
even though the surface of the metal was given 
only an ordinary cleaning. Because of the diffi- 
culty of removing replicas from slightly hydro- 
phobic surfaces, however, it is recommended that 
the cleaning process described previously be 
carried out when making replicas with dioxane 
and polyvinyl formal resin. 

Several methods>** have been described for 
preparing thin films for holding various types of 
specimens for study in the electron microscope or 
electron diffraction camera. 

The following procedure has been developed in 
this laboratory and will be described since it 
permits a somewhat better control of thickness 
and conditioning treatment than the other tech- 
niques known to us: 

A glass plate, such as a microscope slide, 
is cleaned by rubbing with MgO powder or 
with a standard glass cleaning solution such as 
K»,Cr.0;—H2SO,. After drying, the slide is 
dipped into a 0.1-percent solution of polyvinyl 
formal in ethylene dichloride or dioxane, and 
withdrawn and the solvent evaporated. After the 
solvent evaporates (a matter of less than a 
minute), the film is scored with a needle or similar 
sharp pointed tool around the edges of the plate 
and at the same time the main part of the film 
may be subdivided into as many smaller areas as 
may be desired. The plate is then held vertically 


®L. Harris and E. A. Johnson, Rev. Sci. Inst. 4, 454 
(1933). 
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Fic. 6. Sample of colloidal V:0; mounted on conditioned 
surface of polyvinyl formal film. 


over a tray containing water whose surface is 
clean and with a quick motion plunged under 
water. The plate is then withdrawn with a slow 
motion to permit the water to peel from the sur- 
face and again lowered, this time at a somewhat 
slower rate. The film will peel off the plate and 
float on the water surface. 

We have found that freshly split mica, prefer- 
ably muscovite, as a base for forming thin films 
has certain advantages over glass slides since such 
a surface does not require cleaning and is very 
smooth. A fresh splitting should be used each time. 

For mounting the film an ordinary ?” or 3” 
steel washer is given a thin coat of adhesive such 
as rubber cement, held in a horizontal position 
with a small magnet, and lowered into contact 
with the floating film. The mounted film is then 
removed from the water by raising the frame at 
one edge and progressively tipping it upward and 
away from the water. As this is done, the water 
peels off the under surface of the film and it is 
removed completely dry and free of any con- 
taminating residue or internal strains. In this 
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manner, using polyvinyl formal, a film only 100A 
in thickness and an area of 3 cm? is easily ob- 
tained. Such a film reflects little light and is 
similar in appearance to a “‘black’’ soap bubble. 
A number of these thin films can be prepared at 
one time and kept for future use. 

By placing specimens which tend to fly off 
when subjected to the electron beam between 
two films in contact with each other, they are 
securely anchored and remain in place in the 
field. The smoke particle illustrated in Fig. 5 is an 
example of a specimen which is difficult to ob- 
serve unless held in place between two films. 
When specimens in water solution are to be 
mounted, it is sometimes advantageous to have 
the film hydrophilic. This is easily accomplished 
by forming a monolayer! of a protein, such as 
insulin or egg albumin, on a water surface, 
putting it under a pressure of a few dynes cm, 
lowering a mounted film into contact with it, and 
then removing the film again from the water sur- 
face. This process produces a surface which is 
completely wet with water with an increased 
thickness due to the deposited film of not more 
than 30A. A few drops of the liquid containing 
the specimens to be examined are placed in con- 
tact with the wet surface and then the excess 
water is drained and the film dried. Figure 6 
illustrates a sample of colloidal V2O; prepared in 
this manner. 

Replicas may be stored under ordinary atmos- 
pheric conditions for a long time. In our collec- 
tion, we have polyvinyl formal films 100A in 
thickness which are now nearly a year old and 
others less than 1000A which show no sign of 
deterioration after more than three years. Some 
of our steel surface replicas used for demon- 
stration purposes have been in and out of the 
electron microscope more than a score of times 
over a period of several months without showing 
any ill effects. 

In conclusion, it might be pointed out that we 
find it possible with the equipment and technique 
described in this paper to have a good image of 
the surface of a metallographic specimen on the 
fluorescent screen of the electron microscope 
within five minutes after the sample has been 
received from the polishing room. 


10[, Langmuir and V. J. Schaefer, Chem. Rev. 24, 
181 (1939). 
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The Matrix Theory of Torsional Oscillations 


Louis A. Pipes 
Graduate School of Engineering, Harvard University, Cambridge, Massachusetts 
(Received March 5, 1942) 


The general case of vibration of a shaft with several rotating masses is discussed. The 
analysis is reduced to a simple matrix multiplication of second-order matrices and the natural 
frequencies of symmetrical systems are discussed. The connection of the matrix method and 
the well-known Holzer method of computing the natural frequencies of torsional oscillations 
is considered. The analysis of geared systems and continuous shafts is discussed. 


INTRODUCTION 


HE analysis of problems on_ torsional 

vibrations in electrical machinery, Diesel 
engines, and propeller shafts can be reduced to 
a system of several rotating masses on a shaft.” 
In the usual analysis of the torsional oscillations 
of the crankshaft of an internal combustion 
engine, the actual system is replaced by an 
equivalent system consisting of a uniform shaft 
carrying equidistant identical disks. The regu- 
larity of the system leads naturally to an analysis 
using the calculus of finite differences.*~* In this 





discussion the problem will be analyzed by the 
use of simple matrix multiplication. The matrix 
analysis emphasizes the physical significance of 
the mathematical procedure and appears most 
direct, since the problem of the evaluation of 
troublesome arbitrary constants is not en- 
countered. Geared and continuous systems may 
be most effectively analyzed by the matrix 
method. The connection between the well-known 
Holzer method of analysis‘ and the matrix 
method is established. 

1J. P. Den Hartog, Mechanical Vibrations (McGraw- 

Hill, 1940), Chapter 5. 
' 2F, M. Lewis, “Torsional vibration in the Diesel 
engine,”’ Trans. Soc. Nav. Arch. Marine Eng. 33, 109 (1925). 

*M. A. Biot, ‘Equations of finite differences applied to 


torsional oscillations of crankshaft,” J. App. Phys. 11, 
530-537 (1940). 

‘'M. A. Biot, “Vibrations of crankshaft-propeller sys- 
tems, new method of calculation,” J. Ae. Sci. 8, 107-112 
(1940). 


5 Th. von Karman and M. A. Biot, Mathematical Methods 
in Engineering (McGraw-Hill, 1940). 
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I. STATEMENT OF THE PROBLEM 


The problem under consideration in this dis- 
cussion is the determination of the natural fre- 
quencies and modes of oscillation of a set of disks 
on a shaft, as shown in Fig. 1, when the torsional 
elasticity of the shafting c, and the moment of 
inertia I, of the disks are specified. 

The assumption will be made that the disks 
are rigid members and that the connecting 
sections of shafts are devoid of inertia. The ter- 
minal disks will be assumed to be either free or 
fastened to fixed supports through elastic shafting 
in a specified manner. The case of geared systems 
and systems where the inertia of the shafting 
must be taken into account will be considered 
also. 


II. THE ASSOCIATED MATRICES OF DISKS 
AND SHAFTS 


a. Associated Matrices of a Disk 


Consider a disk of moment of inertia J about 
its central axis AB as shown in Fig. 2. Let M, 
be a couple applied to the left side of the disk 
tending to cause the disk to rotate in a clockwise 
sense about the axis AB. Let My, be a couple 
applied to the right side of the disk tending to 
cause the disk to rotate about AB in the counter- 
clockwise sense. Let ¢: be the angular displace- 
ment of a point on the left side of the disk with 
respect to some fixed point in space and let ¢2 
be the angular displacement of a corresponding 
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Fic, 2. Stiff disk. 
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point of the right-hand side of the disk. By 
Newton’s law for a rotating body, we have: 


M,—M.,=I1(d*$,/dt*). (1) 


Since the disk is assumed devoid of elasticity, 
we have: 


o1= 2. (2) 


Now let it be assumed that the couples applied 
to the disk are of an oscillatory nature of angular 
frequency. We then have 


M,= T, sin wl, (3) 
M.= T2 sin wl. (4) 


Let us assume also that these oscillatory 
couples give rise to an angular deflection of an 
oscillatory nature. Hence: 


oi1= 6; sin wl, (5) 
o2= Ae sin wt. (6) 


Substituting these relations into Eqs. (1) and 
(2) and dividing out the common factor sin wt, 
we obtain 


T1=T2—Iu62, (7) 
= A>. (8) 


Equations (7) and (8) may be written in the 
convenient matrix form: 


T | 1-—w°J71| 72 
lello itlat = 
a) to 1Jle 
The associated matrix 
1—w*l 
ra (10) 
0 1 
will be called the associated matrix of the disk. 
b. The Associated Matrix of a Shaft 


Consider a section of shaft as shown in Fig. 3. 


A ----| _ |e--- B 


Bb; ¥, 
Fic. 3. Massless shaft. 








Let the torsional stiffness of the shaft be c. 
Assume that the inertia of the shaft may be 
neglected. Let M, be a couple applied to the 
shaft tending to turn the left end of the shaft in 
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a clockwise sense about the central axis of the 
shaft A—B. Let Mz be a couple applied to the 
right end tending to turn the right end of the 
shaft in a counter-clockwise manner about its 
central axis A—B. Let ¢; be the angular displace- 
ment of a point on the left end of the shaft with 
respect to a fixed reference point in space, and 
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Fic. 4. Combination of disk and shaft. 





let @2 be the angular displacement of a corre- 
sponding point on the right end of the shaft with 
respect to the fixed reference point. Let ¢:=¢2 
when the shaft is not undergoing torsional stress. 
We then have by Newton’s third law: 


M,=M:2. (11) 
By the definition of torsional stiffness, we have: 
c($1— $2) = M2. (12) 


Now, as before, assume that the couples that 
are applied to the shaft are of an oscillatory 
nature and are given by (3) and (4). Let the 
angular displacements be given by (5) and (6). 
Substituting these equations into (11) and (12) 
and dividing out the common factor sin (wt) we 
obtain: 


T1=T;, (13) 
#=(Ts/c)+by (14) 


These equations may be written in the con- 
venient matrix form: 


Ts 1 O72 
lallive at 
0) Life 1) 6 


The matrix 
1 0 
(ul. -|. | (16) 


e 7% 


will be called the associated matrix of the shaft. 


c. Associated Matrix of a Disk Attached to 
a Shaft 


Consider the system of Fig. 4 which consists 
of a stiff disk of moment of inertia J attached to 
a massless shaft of torsional stiffness c. Let J; be 
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the oscillating torque applied to the left side of 
the disk in the clockwise sense and 7, the oscil- 
lating torque applied to the right end of the shaft 
in the counter-clockwise sense. 6; and 6. are the 
reference angles of the left side of the mass and 
the right end of the shaft, respectively. In view 
of (9) and (15) we have: 


| 2] } ~ 1 Ae - 

= z (17) 
las) Lo 1 Jdhi/e 15 0! 

This relation follows from. the fact that the 
torque transmitted through the right side of the 
disk is equal to that impressed on the left end of 
the shaft and the angular displacement of the 
right side of the disk is equal to that of the left 
end of the shaft. 

Multiplying the two square matrices of (17) 
together we obtain: 


Area 


Ill. FREQUENCY EQUATIONS OF SIMPLE 
SYSTEMS 


a. Disk and Shaft Attached to Fixed Support 


(18) 


The associated matrices of simple systems are 
of great value in determining the frequency 
equations of these systems. Let us consider 
first the simplest system of Fig. 5. The frequency 
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equation giving the free oscillations of this 
simple system may be obtained from (18) by 
noting that at the fixed support @2=0, and that 
for free oscillations 7,;=0. In this case we have 


directly from (18) 


T,=T.(1—w*I/c) =0. (19) 
Hence 
(1—w?I/c) =0, (20) 
or 
w=(c/I)}. (21) 


A slightly more complicated system is that of 
Fig. 6. This system consists of two disks of 
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moments of inertia J; and Jz coupled by means 
of two elastic shafts (as shown) to a rigid 
support. The equations have a somewhat simpler 
form if we write the associated matrices of the 
system in terms of the torsional flexibility con- 
stant k=1/c rather than the torsional stiffness c. 
In this case, we have #2:=0 since the shaft is 
fixed at the right end. Consequently we have: 


("| eer 
a) ky 1 


(1—w*leks) —w* ls { 72] 
| | . (22) 
ks 1 jlo) 
Performing the matrix multiplication, we 


obtain: 
T= T>2} (1 —w’*I,k,)(1 — wok) —w* Ik! az @), (23) 


Since the left end is oscillating freely, 7,=0 and 
hence the expression in brackets must vanish. 
This leads directly to the frequency equation: 


wT Tekike— w*(TikitTeke+like) +1=0. (24) 


As another example of this method of obtain- 
ing the frequency equation of the system, let us 
consider the system of Fig. 7. In this case ki, ke, 
k; are the torsional flexibilities of the shafts as 
shown and J, and J, are the moments of imertia 
of the disks about the axis of rotation. 

In this case 6.=0 since the system is fixed at 
the right end. Therefore we may write the 
relation 


| } \ reed 
los Le 1 ke 1 


(1 —w*Tok3) —w*l21{ T2 
x| || (25) 
k3 1 0 


Performing the matrix multiplication, we 
obtain: 
6; = T.[ki(l —w*I,k2)(1 —w*Iok3) 
—w* Ti kik3t+ke(1 —w*Iok3) +k3 |=0. (26) 
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Since the left end of the system is rigidly 
fastened, 6,=0. It follows that the term in 
brackets must vanish. This leads to the frequency 
equation: 


wT Tok ikok3— wi (kiko +hkiks) 
+1o(Rkikst+koks)+(Ritkotks) =(0 (27) 


for the frequency equation of the system of Fig. 7. 

As a final example of this method of determin- 
ing the frequency equation of torsional vibrating 
systems, let us consider the system of Fig. 8. In 
this case both ends of the system consisting of 
three disks and two shafts are free. Consequently 
we have 7:=0. We may, therefore, write the 
relation: 


| 71) (1 —w* ky) — wT) 
lea! Lose, 1 








(1 — wT oko) —w* I) 1 —w°' ls | 0 
| | f (28) 
.” * 1 ilo 1 dle 


Performing the matrix multiplication, we 
obtain: 


7 0.= —w*I;[ (1 —w*Tik,)(1 —w* Tsk) —w Tike | 
—w'ls(1—w*I1k1) —w*T;. (29) 


Since the left end of the system is free, 7,=0, 
so the right member of (29). must vanish. This 
leads to the conditions: 


w= 0, (30) 
or 


wT ToT 3kike—w*[ Ti Tokit DiI sk 
+ TiIsko+TeIsk2 J+ (it+12+Is) =(. (31) 


More complicated systems having a greater 
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number of disks may be analyzed in a similar 
manner. 


IV. FREQUENCY EQUATION OF EQUIDISTANT 
IDENTICAL DISKS 


In the analysis of the torsional oscillations of 
internal combustion engines, the vibrating 
system is reduced to one of a uniform shaft 
carrying equidistant identical disks. The analysis 
of systems of this type has been carried out by 
Biot*-* using the calculus of finite differences. In 
this section the analysis of these symmetrical 
systems will be considered by the use of matrix 
algebra. 

To fix the ideas, let us consider the oscillations 
of the system of Fig. 9. This is a system of n 


Fic. 9. 


identical disks of moments of inertia J about the 
axis of rotation fastened together by sections of 
shafting of torsional flexibility k and free at both 
ends. The system of Fig. 9 may be regarded as 
made up by the interconnection of m fundamental 
units of the type shown in Fig. 10. 


I 
/2 a 


= 














Fic. 10. Fundamental unit. 


This fundamental unit consists of a disk of 
moment of inertia J to which are fastened two 
sections of shafting of torsional flexibility k/2 as 
shown. The associated matrix of this fundamental 
unit is given by: 


[ Lo. ‘| 
k/2 ifO 1 SLR/2 1 


(32) 








Now since the system of Fig. 9 is a chain of n 
units of the fundamental type of Fig. 10, we may 


write: 


Ti A Bry{T>» 
fe Tos 
6; C D A. 
where A, B, C, and D are the elements of the 


associated matrix of the fundamental unit given 
by (32). We note that 


A=D=(1—o°Ik/2). (34) 


This is a consequence of the symmetry of the 
fundamental unit. 

To obtain the frequency equation of the 
system of Fig. 9, we could, of course, raise the 
associated matrix of the fundamental unit to a 
high power and then follow the procedure of 
Section III. This would be a very tedious pro- 
cedure. Instead we take advantage of a theorem 
in matrix algebra which enables us to write down 
the power of a matrix directly. 

This theorem is established in the Appendix. 
The theorem states that 


cosh (an) Zo sinh (an) 


A By" 
= | sinh (an) , (35) 
C D — cosh (an) 
0 
where 
A+D w*Ik 
cosh (a) =— =(1- ). (36) 
2 2 
((A+D)?—4)! 
es 2c . 
sinh (a) sinh (a) 
= ——_=——_——,_ (37) 
k cosh? (a/2) c 
provided that 
A=D (38) 
and 
(A+D)+#¥ +2. (39) 


To simplify the writing let us use the notation 


cosh (an) Zo sinh (an) 


sinh (an) 
———  _ cosh (an) 
Zo 


A, B, A By 
-| |-| ]- a 
C, D, Cc D2 
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a. Natural Frequency of Free Shaft 


We are now in a position to compute the 
natural frequencies of the free shaft of Fig. 9. 
Since the right end of the system is free, 7;=0 
and we may write: 


("| * aie 
A; 7 Ge D,, Ao , 


Carrying out the matrix multiplication, we 
obtain: 


(41) 


7, =B,,02. (42) 


Since the left end of the system is free, 7,=0 
and it follows that: 
B,,=Z, sinh (an) 

=[sinh (a) sinh (an)/k cosh? (2/2) ]=0. (43) 


This is the frequency equation of the system 
of Fig. 9. From this we obtain 


sinh (an) =0, (44) 
and 
an=rrj, r=0,1,2,--- (45) 
Hence 
a,=rxj/n, r=0,1,2,°--, j=vV—-1. (46) 


Now the condition (39) excludes the value 0, 
and n that may be taken by rv. The value r=0 
corresponds to the system rotating as a whole 
without any relative motion with the system. 
To obtain the natural frequencies we have 


Ik 
cosh (a,) = ( -u—) 
2 
rm] 1K 
= cosh ( - ) =0os ( ). (47) 
n n 
Hence 
2 ‘Tr 4 as 
wai — cos ()|-- — sin? ( "), (48) 
Ik n Ik 2n 


or 
2 ‘Tr 
-—_ sin (~*) r=1,2,-+-(n—1). (49) 


The possible natural angular frequencies of the 
system are given by (49). 
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b. Natural Frequency of Shaft Fixed at 
Both Ends 


Let us now consider the system of Fig. 11. In 
this case the two ends of the system are fastened 


1 2 3 T. 
6, : «&£ -« I 62 
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to rigid supports as shown. We can write #.=0, 
and hence 


(7 1 OyA, Bayi O77? 
CO Pa A 
load Le2 idle, D,Jjle/2 1Jlo 
Carrying out the matrix multiplication, we 
obtain: 


6:/T2=(Ant+D,)k/2+Ci+Bak?/4. (51) 


Since the left end of the system is fixed, #:=0 
and hence the right member of (51) must vanish. 
Substituting the values of A,, Bn, Cr, and D, 
from (35) into the resulting expression, we 
obtain: 


cosh (an)+sinh (a) sinh (an) 


4 cosh? (a/2) 


sinh (an) cosh? (a/2) 
+—____—_———-=0, (52) 
sinh (a) 
or 


[sinh (n+1)a/sinh (a) ]|=0. (53) 
The solution of the transcendental equation 
(53) is 
a,=rmrj/(n+1), r=1,2,3,---n. (54) 
Substituting these values of a into (36), we obtain 
2 Tr 


w= — «in| "|. r=1,2,-+-n (55) 
(Ik)} 2(n+1) 


for the natural frequencies of the system of Fig. 
11. 


c. Natural Frequencies of Shaft Fixed 
at One End 


The system of Fig. 12 represents a uniform 
shaft carrying identical disks equidistantly 
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6, 
Fic. 12. 


spaced, free at the left end and fixed at the right 
as shown. In this case we have 62=0 and we can 


write: 
[71 4A, Bz 1 OV 72 
Lalle olla allot 
6; ce. D,, k/2 1 0 


Carrying out the matrix multiplication, we have: 
k 

Ti= (4.48. )T (57) 
Since 7,=0, we must have 


k 
(4 +B.) =cosh (an) 


sinh (a) sinh (an) 








—-- =(0. (58) 
2 cosh? (a/2) 
This may be written in the form 
(2n+1)a a 
cosh — cosh (<) 
2 2 
ed eA 
cosh? (a/2) 
Hence 
cosh [(2n+1)a/2]=0. (60) 
This leads to 
(2r—1)xj 
a, =———__,, r=1,2,---n. (61) 
(2n+1) 


Then from (36) we obtain 


2 (2r—1)x 
w,= -sin| —| y=1,2,---n. (62) 
(Ik)? 2(2n+1) 





d. Natural Frequencies of General 
Terminal Conditions 


In the simple cases that we have considered, 
the frequency equations have been particularly 
simple transcendental ‘equations whose solution 
was easily written down. Let us now consider 
the case of Fig. 13. This arrangement may be 
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considered formed from the system of Fig. 9 by 
the addition of shafting of torsional flexibility. 


k,=K,-—k 2, 
ko=K.—k y a 


(63) 
(64) 


at the left and right ends, respectively, and two 
disks of moments of inertia J; and Js as shown. 





rn t 
1 Iz 2 
[ROL L -4 
1 
8, 2 3 n 6, 
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By the methods of Section III we can obtain the 
associated matrices of the terminal portions of 
the system by direct matrix multiplication. We 
may then write: 


(71) ewes 
| a! ky 1 


A n B,, 1 —_ wT» Ts) 
le. ol Ita 
Ce D, ke (1 —w*I ok) 6.) 


(65) 


For free oscillations of the system we must 
have 7,;=0, 7:=0. Carrying out the matrix 
multiplication, we obtain after some reductions 


tan (bn) =[M,M/(M2M?— M3) ], (66) 
where 
a=(jb), j=V-1 (67) 
and 
sin (b) 
M= —, (68) 
k cos? (b/2) 
M,=o0'T Jo(ki tke) —w*(114+- 12), (69) 
Me=wo'l Tok ike—w?(11ki +-Tok2) +1, (70) 
M;=0'l le, (71) 
2 b 
w= sin (-). (72) 
(Ik)} 2 


By means of (72), the right member of (66) is 
transformed to a function of 6. The tran- 
scendental Eq. (66) may be solved graphically 
for 6. From these values of b and Eq. (72) we 
obtain the natural frequencies of the system. A 
form of Eq. (66) more suitable for graphical 
solution has been obtained by Biot.’ 
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e. Connection with Holzer’s Method 


H. Holzer® has given a method for the calcu- 
lation of the natural frequencies of torsional 
oscillations of shafts. In the interest of unity it is 
well to show the relation of the matrix method 
to Holzer’s method. To fix the ideas, a system 
discussed in reference 1, page 227, will be con- 
sidered. This system is shown in Fig. 14. 





Fic. 14. 


This system represents the equivalent of a 
six-cylinder Diesel engine with a flywheel directly 
coupled to a ship propeller through a long pro- 
peller shaft as regards torsional oscillations. It is 
evident that this system may be regarded as 
formed by the interconnection of six identical 
fundamental units, as shown in Fig. 10, a section 
of shaft of torsional flexibility k/2, a disk of 
moment of inertia J;, and a system consisting of 
a shaft of torsional flexibility k:, and a disk of 
moment of inertia J, Accordingly, we. may 
write: 

* b “| 1 | =o 
lat Le, Dedle/2 ttlo 1 - 


1 O1f 1—w"l2 [22 
x| ; (73) 
k, itlo 1 dle, 


out the matrix multiplication, we 





Carrying 
have: 


T, As Bs 
| | -| k w*Iik 
2 2 
[= wT» T > 
x| | ; (74) 
ky ( 1 —w*lI,k, A> 
By introducing the transformation (67) into 
Eq. (40) we obtain: 


wl, 


A,=cos (bn) =D,, (75) 


®H. Holzer, Die Berechnung der Drehschwingungen 
(J. Springer, Berlin). 
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—sin (bn) sin (bd) 


2 erry, (76) 
k cos? (b/2) 
sin (bn) b 
C,=———_k cos? (-). (77) 
sin (b) 2 


From (72) we have the relation 


b w(Ik)? 
sin (-)- —~, (78) 
2 2 


The essential feature of the Holzer procedure 
is to realize that in the case of free oscillations 
T,=T.=0. If then, one assumes an arbitrary dis- 
placement 62 having placed 7.=0, then if w in 
(74) has the value of one of the natural fre- 
quencies, 7; will come out equal to zero. The 
method of Holzer consists in trying various 
values of w in Eq. (74) until one is found that 
gives 7,;=0 for an arbitrary choice of 62, rather 
than in solving the transcendental equation 
which results when (74) is expanded and 7; set 
equal to zero. 

The advantage of formulating the Holzer 
process in matrix form is that if the system 
possesses component parts of identical sections, 
then the associated matrices of these parts may 
be written down and the numerical work greatly 
reduced. For example, in the system of Fig. 14 
it is evident that if the system is vibrating in the 
fundamental mode, the system may be replaced 
by a system of two disks connected through a 
shaft for an approximate calculation of w; the 
fundamental angular frequency. Starting with 
this first value of w we use (78) to calculate 0 
and (75), (76), and (77) to compute Ag, Be, Ce, 
D,. Then we see whether we obtain 7,;=0 when 
T.=0 and @, is arbitrary. By successive trials we 
adjust w until we do obtain the required relation; 
this gives w,, etc. In (3) the natural frequencies 
of the system of Fig. 14 are obtained by using 
the methods of finite differences. 


V. GEARED SYSTEMS 


The matrix method is very well suited to the 
analysis of torsional oscillations of systems whose 
component parts contain gears and pinions. All 
that is necessary is to obtain the associated 
matrix of a system composed of two massless 
gears, as shown in Fig. 15. Let the gear ratio of 
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1 
Fic. 15. 


the upper gear to the lower gear be n: 1. We 
then have the relations: 
6.= —nh,, T2.= = (T;/n), (79) 


or in matrix form 


T; = 0 T2 
Sieca:2 
6, 0 —(1/n)Jl 0 


The square matrix of Eq. (80) is the associated 
matrix of the gear system of Fig. 15. As an 


T, I, 


Ry 














6, 
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example of the use of this associated matrix, let 
us obtain the frequency equation of the system 
of Fig. 16. Let the gear ratio of the upper gear 
to the lower gear be m and let the moments of 
inertia of the two gears be 7; and 72 as shown. The 
two gears are connected to disks of moments of 
inertia J, and J, by shafts of torsional flexibility 
k, and k2 as shown. We then have: 


"| eater. 
a) bk 1 


1 —w"t; [_—n 0 1 —w'ts 
4, he -aniies! 
0 1J4L0 -(1/n)ilo 1 


v1 _ wl T>2 
x || | (81) 
ike = (1—w*lok2)d! A. 





Setting 7,.=0 and performing the indicated 
matrix multiplication, we obtain from the equa- 
tion 7, /6.=0, the following frequency equation: 


w' Tol Tok ike —w*To(I:ki + Toke) 
—w*(I,Ton?k,— I, Tok2) + 10+ (n?Io+J;) =0, (82) 
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where 


To = (n*i2+1;). (83) 


More complicated systems may be discussed in 
a similar manner. 


VI. CONTINUOUS SHAFTS 


In the discussion of the previous sections, the 
mass of the shafts of the systems considered was 
neglected. Let us now consider the case where the 
mass of the shaft, and hence its moment of 
inertia, may not be neglected. In such a case, we 
may regard the shaft as being made up of an 
infinite number of infinitesimal sections of the 
type of Fig. 17. 


I Ox 





Ros a 
<—_ — > 


Ax 
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The associated matrix of this fundamental 
unit is given by (32) and may be written in the 
form: 


r w* TkAx* ; 
(:- ) _ w*TAx 
2 
k*w*TAx' w*TRAx? 
pny 6-2) 
4 4 2 d 


A B 
= , (84) 

Cc Dp 
where J is the moment of inertia per unit length, 
and k the torsional flexibility per unit length of 

the shaft. 

The procedure is to pass to the limit Ax—-0 
and consider the continuous shaft as composed 


of an infinite number of sections of the type of 
Fig. 17. From (36) we have 








w*TRAx? 


A+D 
2 2 


cosh (a) = 


a‘ 


a’ 
=1+—+—+-:--. (85) 
2 4! 
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From this we obtain the relation 


lim a=jw(Ik) Ax, 


Ax “0 


where j= \V—-1. (86) 


Now if is the number of sections of the type 
of Fig. 17, we have: 


lim n-Ax=1, 
n—-2 
Ax—0 


(87) 


where / is the length of the shaft. Now from (37) 
we have 


w*TkAx*\? ; 
it->) 4 
2 


lim Z= 
Ax—0° 2kAx (88) 
. I ; . 
= iol ) j=yv-1. 
k 
If we now introduce the notation: 
¢=wl(Ik)', (89) 


then by making use of (86), (87), and (88) and 
substituting into (35), we obtain the associated 
matrix of a shaft of length /, moment of inertia 
per unit length J, and torsional flexibility per 
unit length k to be: 


cos @ —ysin @ A B 
| lca * 
sing/y cos¢@ C D 


where 


(91) 


As an example of the use of the associated 
matrix of a uniform shaft, let us calculate the 
natural frequencies of the system of Fig. 18. The 


T, L T, 

















1 62 
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system of Fig. 18 consists of a uniform shaft with 
two disks attached to its end. We are interested 
in the free vibrations of the system. In this case, 
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we may write: ; 
| “| : = cos @ —y sin *] 
| 0; 0 1 sin ¢/y cos@ 


1—w’l, 0 
| | ; (92) 
0 1 Os 


Carrying out the matrix multiplication, we 
obtain: 


T, 
—_= —w?l, cos o-y sin o 
62 
w*I,I, sin 
——w’*I,cos¢=0. (93) 





y 
This reduces after some reductions to 


wy(I,+J2) 
tan ¢=———_—__. (94) 
wl saat y? 


If we now let Jy)=lI be the total moment of 
inertia of the shaft, and introduce the notation: 


m=I1;/Io, n=I>/Ipo, (95) 


Eq. (94) may be written in the form: 
tan ¢=—————__. (96) 


From (89) we obtain the relation: 
w= o/l(Ik)}. (97) 


Equation (96) is a transcendental equation 
for @ and the various angular frequencies are 
obtained from the roots of this equation by (97). 
The solution of this equation may be easily 
carried out graphically for any special case. If 
the left end of the system is fixed, then J; < 
and (96) reduces to 


tan ¢=1/nd¢. (98) 


If the moments of inertia of the disks J, and J, 
are much larger than Jo, then (96) reduces to 


@ tan ¢=(1/m+1/n). (99) 
Many other special cases may be considered. 
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VII. CONCLUSION 


The simple and natural manner in which 
matrix algebra enters into the analysis of tor- 
sional oscillations is evident from this discussion. 
The conciseness of the results and the unity of 
approach should prove helpful to workers in this 
field. The simple operation of matrix multi- 
plication replaces the more complex methods of 
differential equations, Holzer tables, and dif- 
ference equations in the conventional treatment 
of the subject. The only advanced result used is 
Eq. (35). This is derived in the Appendix. 


APPENDIX 


Consider the matrix 
(1 -) a; 
2 A B 
(uJ -| | (1 
(s-“2) (: _) Cc 2 
4 2 


The determinant of the matrix [u] has the important 
property that 





A Bi 
lu|= | =AD-BC=1. (2) 
C BD 
The characteristic matrix of [u] is given by 
(q-A) -—B 
Cf(g)J= : (3) 
-C (q-D) 
and the characteristic equation of [u] is 
A(q) = | f(g) | =g—q(A+D)+1=0. (4) 
The latent roots of [u] are, therefore, 
A+D A+D)? ; 
o.2=(-)+[S*-1] , (S) 
provided that 
(A+D) ¥ +2. (6) 
The latent roots have the property that 
qiq=l, | (7) 
gat+gq=(At+D). (8) 


In the special case we are considering, the matrix [u] 
has the property that A =D. Now let 


qi=e«*, g2=e%. (9) 
Then in view of (8) we obtain 
cosh(a) =(A+D)/2=A. (10) 


The adjoint matrix of [u] is given by 


(q—D) B 
[F(q)]= . (11) 
Cc (q-A) 
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Now by Sylvester’s theorem’ we have that provided q: 
and gq, are distinct 


u\= 
2qi—(A+D) 84 (qi—A) 


qo" (q2—D) B 
+ P . (12) 
2q2—(A+D) rt (q2—A) 


Now on making use of (2), (10), and the fact that 


™R. A. Frazer, W. J. Duncan, and A. R. Collar, Ele- 
mentary Matrices (Cambridge, 1938), pp. 78-79. 


(A =D) we obtain after some reductions, the relation 


cosh (an) Zposinh (an) 


(u}"=| sinh (an) — (13) 
—————— cos tf 
Zo 
where n is a positive or negative integer, and 
Zo=sinh (a)/C. (14) 


This is Eq. (35) of Section IV. Applications of this 
theorem to electrical circuit theory are found in reference 8. 

8L. A. Pipes, ‘“‘The matrix theory of four-terminal 
networks,’’ Phil. Mag. 30, 370 (1940). 





The Mechanism of the Electric Spark 
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Meek’s criterion for the formation of a spark discharge 
in a gas is criticized. An error in the derivation of his 
formula for the radial field about an electron avalanche 
channel is noted, and the value adopted for this field is 
shown to be too small to draw the postulated number of 
electrons into the channel from the space outside. The 
importance attached to the radial field is found to be 
misplaced because of the above and because a criterion 
based solely on the total number of ions produced in an 
avalanche gives results for sparking potentials for different 
spark lengths and for a gas at different pressures in 
agreement with those obtained by Meek’s method. Both 
methods give results which, using the values assigned to 
the constants, fall more and more below the experimental 
values of the sparking potentials as the length of the 


INTRODUCTION 


HE application of the Wilson cloud chamber 
to the study of the distribution of ions 
along the paths of disruptive discharges at 
different, very short time intervals after the 
discharge voltage is applied to the gap revealed 
some interesting facts about the early stages of 
such discharges. In discharges from pointed 
conductors to which such study was first applied,' 
the ion tracks for negative discharges always 
show the ions distributed throughout a single 
conical-shaped volume rounded at its forward 
end, whereas with positively charged points the 
tracks, which are more difficult to obtain, often 
show a number of branches. 
tL, B. Snoddy and C. D. Bradley, Phys. Rev. 45, 432, 


751 (1934); 47, 541 (1935); H. Raether, Zeits. f. Physik 
94, 567 (1935). 
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spark gap, 4, is increased until at 20 cm the deviation is 
over 7 percent. To bring the computed value for 6=20 cm 
into agreement with the experimental value the constants 
in the equations would have to be increased by a factor 
210'* over the values which give agreement between 
formulas and experiment for 6=1 cm. The dominant 
exponential term in each formula is based on the incorrect 
assumption that the field is not affected by the space 
charges in the avalanche. That a large amount of ioniza- 
tion, postulated by Loeb and Meek, is produced by a 
high potential wave which travels through the channel as 
soon as the positive streamer reaches the cathode is 
questioned because the fields found in such waves are too 
small to produce ionization by electron collision in air at 
atmospheric pressure. 


In uniform fields the tracks for the shortest 
times after the voltage is applied show a rapid 
development from the cathode of Townsend 
electron avalanches, and these tracks have in 
general the same appearance as those observed 
with negative point discharges. Contrary to 
what would result from electron diffusion alone, 
some of these avalanches cease to increase in 
diameter as they proceed, and indeed appear to 
contract during the last half of their length.? 

The tracks formed in uniform fields after the - 
electron avalanche has reached the anode indi- 
cate that a form of discharge is now moving 
from the anode to the cathode, and these tracks 
possess branches similar to those found on 
tracks with positive discharges from points. 


2H. Raether, Zeits. f. Physik 37, 560 (1936); 107, 91 
(1937). 
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When several electron avalanches are de- 
veloped at the same time the evidence shows 
that the final spark path is evolved from but 
one of these avalanches.? 

Several hypotheses have in recent years been 
advanced to account for the transformation of 
an electron avalanche into a continuous con- 
ducting path between the electrodes. In one of 
these, thermal ionization arising from heat 
generated in the avalanche was invoked. In 
another, the applied field was supposed to be 
augmented so greatly near the cathode owing to 
space charges as to make ionization by positive 
ions possible even in gases at atmospheric 
pressure. Such space charge fields have been 
called upon, too, to produce the essential new 
supply of electrons by detaching these from 
ions which were formed in the gas before the 
avalanche started, or by having high speed 
positive ions eject them from the surface of 
the cathode. 

At the present time the prime agent for 
supplying the necessary electrons to the space 
behind the head of the avalanche is generally 
considered to be photo-ionization. 

The radiation responsible for this ionization 
is believed to have its origin in the avalanche 
itself where some of the molecules have been 
raised to excited states by the high speed 
electrons. That a high frequency radiation 
capable of ionizing a gas is emitted along a part 
of the path of a disruptive discharge has long 
been known,’ but only recently has further study 
directed attention to this radiation,‘ and while 
at first it was hesitatingly looked upon as only a 
possible source of the electrons in question it is 
now considered the main source of such electrons. 

The space charge arising from the compara- 
tively immobile positive ions in the avalanche 
channel which previously was invoked for 
activating other processes is now called upon to 
increase the ionizing power of the electrons 
released there by photo-ionization; and it is 
postulated that when this production of ions is 
sufficiently great a retrograde streamer may 
form which, moving toward the cathode, creates 





3 E. Wiedemann, Zeits. f. Elektrochemie 2, 159 (1895); 
J. J. Thomson, Proc. Camb. Phil. Soc. 10, 74 (1899). 

4A. M. Cravath, Phys. Rev. 47, 254 (1935); L. B. Loeb, 
Rev. Mod. Phys. 8, 267 (1936); H. Raether, Physik. Zeits. 
37, 560 (1936). 


VOLUME 13, JULY, 1942 


a conducting path between the two electrodes 
along which the main discharge eventually 
passes. 

The detailed action of such a positive streamer 
was described by Simpson’ in connection with 
the discharge of positively charged clouds to 
earth. However, Simpson assumed that the 
ionization at the head of the positive streamer 
arose from electrons which are normally liberated 
in air by various agencies, and Cravath and 
Loeb® also at first assumed these agencies to be 
the sources of the electrons, when they applied 
the positive streamer theory to ordinary sparks. 
The number of electrons liberated by these 
omnipresent agencies is now considered inade- 
quate for the purpose and their function is 
assigned to the ionizing radiation developed in 
the avalanche itself. 


MEEK’S CRITERION FOR SPARK FORMATION 


Meek’ has postulated a criterion for determin- 
ing quantitatively the condition under which a 
retrograde streamer will form from an electron 
avalanche and develop into a spark. This 
criterion states that such a streamer will form 
when the radial field X, about the positive space 
charge in an electron avalanche attains a value 
of the order of the applied field X,. The supposed 
function of this radial field is to draw into the 
main avalanche channel some of the electrons 
which are liberated outside of it by photo- 
ionization and thus make of the channel a 
conducting plasma. 

Meek evaluated, roughly, the radial field near 
the head of the avalanche in terms of the total 
number of positive ions in the whole avalanche, 
and by reason of the postulated criterion set 
this equal to KX,. Then assuming arbitrarily a 
value unity for the constant K, the field X, 
necessary for the production of a spark was 
found by a trial and error method by means of 
the experimentally obtained relation between 
Townsend's ionization coefficient a and X/p, 
p being the gas pressure. Later the value 0.1 
was .adopted for K because the use of this 
number gave an exact agreement between the 
computed and experimental values of X, for a 





5G. C. Simpson, Proc. Roy. Soc. 111, 56 (1926). 
6 A. M. Cravath and L. B. Loeb, Physics 6, 125 (1935). 
7]. M. Meek, Phys. Rev. 57, 722 (1940). 
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spark gap 6=1 cm, although the computed result 
was not altered much by this large change in K. 

I propose to show that the derivation of the 
formula for the radial field as given by Meek is 
erroneous; that when K is given the adopted 
value 0.1 the radial field cannot be expected to 
draw any large number of electrons into the 
avalanche from the outside, as is postulated; 
and that it is highly illusory to find support for 
the role played by radial fields in spark phe- 
nomena in the fair agreement that exists between 
sparking potentials computed by Meek’s method 
for sparks of different lengths and the corre- 
sponding values found by experiment. 

Some general features of an electron avalanche 
will first be considered. After an electron has 
traversed a distance x in a uniform field the 
number of electrons in the avalanche produced 
by it will be »=e**, where a is the ionization 
coefficient. As actually determined experimen- 
tally m must include not only the electrons that 
owe their liberation by collision directly to the 
initial electron but also those produced anywhere 
in the inter-electrode space by radiation, as well 
as the progeny of such electrons. At the present 
time these additional electrons are thought to 
comprise but a small part of the total number. 
the field remains essentially 
unchanged while an avalanche is advancing, all 
of the electrons in it are confined to a more or 
less spherical volume at its forward end, the 
diameter of which increases with time owing to 
diffusion and the action of electrical forces. The 
density of the electrons in this volume is largest 
at the center of the sphere and diminishes 
outward toward its indefinite periphery. In the 
forward part of the sphere electrons predominate 
greatly over positive ions while at its rear end 
the positive ions are largely in excess. 

‘Raether* found the radius of the head of an 
avalanche in air at atmospheric pressure, 1 cm 
from its origin, to be 0.013 cm when the applied 
field was 30 kv/cm. A calculation shows that 
at that distance, in the volume occupied by 
them, the electrons predominate over positive 
ions in the ratio of about 5 to 1. The remainder 
of the positive ions are in the trail of the ava- 
lanche extending back to the cathode. 


So long as 


* H. Raether, Zeits. f. Physik 107, 91 (1937). 
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The preponderance of electrons in the forward 
portion of the avalanche not only tends to spread 
out the electrons fanwise but also enhances the 
applied field in this region. Since the electrons 
in the avalanche increase rapidly in number as 
it approaches the anode, this space charge 
effect may cause a substantial increase in the 
field near the anode; and as the coefficient a 
increases more rapidly than does the field the 
rate of production of new electrons must here be 
much greater than it is in the unmodified field. 
Such a large increase in the axial field at the 
head of an avalanche may well be the cause of 
what Flegler and Raether® call an instability in 
the avalanches they observed in overvolted 
gaps where a narrower and faster moving 
streamer suddenly shot out from the forward 
end of the avalanche as it approached the anode. 
This was followed by another fast streamer that 
went rearward to the cathode, and which must 
have resulted from an increase in the field in 
this region also, as is generally postulated for all 
positive streamers. That an increase in the field 
owing to the positive charge can occur here is 
conditioned by the fact that the charge is limited 
to a narrow channel and that in this channel the 
charge decreases rapidly toward the cathode. 

Since the fields in front of the avalanche head 
and back of it are both greater than normal 
there must be a region in between the two 
opposite charges where the field is smaller than 
the normal applied field. As the avalanche moves 
toward the anode, these modifications in the 
field should rapidly distort the shape of the 
volume occupied by the electrons from the 
simple spherical form it would have in a uniform 
field. 

Eventually all the electrons in the avalanche 
proper reach the anode, and it is at this final 
stage, when the comparatively immobile positive 
ions alone remain in the channel, that Meek 
proceeds to find the radial field near the anode 
arising from all of the ions in the whole ava- 
lanche. For ease of calculation these ions are all 
assumed to be contained in a spherical volume 
whose radius, 7, is the average distance that 
electrons in the avalanche had diffused from its 
axis in the time taken by the avalanche to cross 


°E. Flegler and H. Raether, Zeits. f. Physik 112, 464 
(1939). 
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the spark gap; and the field X,; is computed for 
this radial distance r. Naturally, many of the 
positive ions were formed at a greater distance 
from the axis than r and these contribute 
nothing to the field at 7. Assuming the ions in 
the spherical volume to be distributed in shells 
of uniform density, the field at a distance 7 from 
the center is, 





ay ’ (1) 


where ¢ is the charge on each ion and 6 is the 
length of the spark gap. 

Meek now introduces the ion density N into 
the formula, obtaining 


X1=(4/3)merN, (2) 


and then eliminates N by a faulty procedure. 
Assuming that the avalanche contains »=e* 
electrons, which presupposes that it has not yet 
reached the anode, he lets each electron advance 
a distance dx so that dn=ae**dx new electrons 
are liberated and an equal number of positive 
ions formed. These oppositely charged particles 
are both distributed over the whole volume of 
the sphere and therefore add nothing to the 
radial field. However, Meek divides this dn, the 
increase in the number of positive ions, by an 
imaginary volume ar*dx, which these ions do 
not occupy, and obtains for the positive ion 
density, N=ae**/rr*, which substituted in Eq. 
(2) gives, 


X 1 = 4ace**/3r. (3) 


This value of X, only differs from the one 
given by Eq. (1) by a factor 4ar/3, since the 
errors in the derivation happen in large part to 
cancel one another. Using values of a@ and r 
given by Loeb and Meek! for the sparking 
distances 0.5 cm and 20 cm, the values of the 
above factor are, respectively, 0.4 and 0.09, and 
these have only a small effect upon the sparking 
potentials computed by Meek’s method. 

Meek attaches great importance to the radial 
field because as already stated he assigns to it 
the function of drawing some of the electrons, 
that are liberated outside of the avalanche by 
photo-ionization, into the avalanche channel and 


electric spark (Stanford University Press, 1941). 
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making of it a conducting plasma. However, 
when K is given the adopted value 0.1, the 
inward deviation of the electron paths produced 
by the radial field is only 0.1 radian at the 
distance r from the axis of the avalanche for 
which X, was derived above, a point which is 
well within the avalanche. For a place somewhat 
farther from the axis, say at a distance 3r, the 
radial field is so weak that electron paths would 
be deviated by it only about a half degree of arc. 
This is certainly without significance for drawing 
enough electrons into the avalanche to make of 
it a conducting plasma, especially as the radial 
field was computed for the part of the avalanche 
near the anode where the electrons have but a 
very short distance to travel before they reach 
the anode. 

It thus appears that radial fields arising from 
positive space charges in the avalanche can play 
only a very minor role, if any, in the initiation 
of a spark. More evidence on this point, quite 
different in character, will be presented later. 

The contention that the radial field is de- 
pendent upon the ion density in the avalanche 
rather than upon their total number m is not 
altogether justified because the appearance of r 
as well as of N in Eq. (2) shows that is also 
involved in X. 


CRITERION BASED ON n=e* 


Now it is the axial fields due to the space 
charges in the avalanche rather than the radial 
fields which must play an important role in 
spark formation since, as already shown, they 
enhance the applied field both in front of and in 
back of the electron cloud in the avalanche. In 
the high applied fields here prevailing the 
ionization coefficient a@ increases rapidly with 
field intensity and hence even a moderate 
increase in the field will increase enormously 
the number of new ions produced by electron 
collisions. 

The only electrons present to the rear of the 
main electron cloud in the avalanche upon which 
the enhanced axial field can act are those 
liberated there by photo-ionization. The in- 
tensity of the radiation which produces this 
ionization is doubtless proportional, at least 
roughly, to the number of ionizing collisions 
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that electrons make in the whole avalanche. 
Since the positive space charge in the avalanche 
is also dependent upon the total number of ions 
in the avalanche, a meaning may thus be 
attached to the fair degree of correspondence 
that is obtained between computed and experi- 
mental values of sparking potentials when it is 
assumed, simply, that a positive streamer can 
establish a conducting path between the elec- 
trodes as soon as the applied field reaches a 
value such as to produce some given large 
number of ions in a single electron avalanche 
while this is traversing the spark gap. 

The simple relation to be now used as a 
criterion for sparking is n=e*’, which assumes, 
as was done for Eq. (1), that the applied field 
remains unchanged by space charges. In evalu- 
ating some sparking potentials by aid of this 
relation the constant value adopted for n was 
10° since this number gives a value of a for 
5=1 cm from which a sparking field in close 
agreement with experiment is obtained, by 
making use of the experimental relation between 
a/p and X/p. The values of a for spark gaps of 
different length are then found from the con- 
stancy of the product aé. 

Calculations of the spark distances obtainable 
in given uniform fields were made by Schumann" 
by empirical equations also based on the con- 
stancy of the product aé, and Loeb and Meek,” 
too, were aware that a formula of this form 
gives results in fair agreement with experiment, 
but they considered the formula as purely 
empirical and so without special significance. 


PASCHEN’S LAW 


Paschen’s law regarding sparking potentials, 
which states that for a given gas the sparking 
potential for a uniform field remains constant 
when the product of the gas pressure p by the 
spark gap length 6 is constant, follows at once 
from the assumption that the liberation of a 
given number of electrons in an avalanche is a 
sufficient criterion for the initiation of a spark. 
Since n=e**=e*?*/”, it follows if m and pé 
remain constant that a/p is also constant and 
hence also X/p and V=(X/p)- pé. 


"W. O. Schumann, Electrische Durchbruchstirke von 
Gasen (Julius Springer, 1923). 
2 Reference 10, p. 118. 
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TABLE I. Variation of sparking potential with air pressure 
p, for 6=1 cm. 


Pp mm Vi kv Va kv 

7600 248.5 248 

3800 131.5 131 
760 32.5 32.2 
380 18.9 18.6 
200 11.35 11.1 
100 6.83 6.55 


Vi—computed for nm =10°. Vy¢—computed by Meek's formula. 

Some sparking potentials obtained in the 
manner indicated for different values of p6 are 
given in Table I under Vj, ” being taken as 10°, 
5 as 1 cm, and p having the values tabulated. 
Except for =100 mm, these potentials are in 
very close agreement with those obtained by aid 
of Meek’s formula," which are shown under V y. 
No experimental data covering this range of 
pressures are available for comparison. 


VARIATION OF SPARKING POTENTIAL WITH 
SPARK LENGTH 


Sparking potentials computed by the two 
methods for air under a constant pressure of 
760 mm when the distance 6 between the 
electrodes is varied from 0.5 cm to 20 cm are 
given in Table II and these are again seen to be 
in good agreement. 

The values under V; are those found when n 
is assumed to be 10° throughout, while under Vy 
are those given by Meek’s formula." 

These computed values are to be compared 
with the experimental potentials, V.x,, found by 
Edwards and Smee!’ for sparks between spheres 
whose diameters were large relative to the 
sparking distances used. 

It is to be noted that, while the computed 
values are in close agreement with the experi- 
mental values for 6 around 1 cm, none the less 
they are farther and farther below the latter 
values as 6 increases, until at 6=20 cm the 
deviation is over 7 percent. The value of 6 for 
which the theoretical and experimental results 
would agree is evidently dependent upon the 
values of K and m assumed in the two formulas 
used. 


18 Reference 10, p. 116. 

4 Reference 10, p. 114. 

16 F, S. Edwards and J. F. Smee, J. A. I. E. E. 82, 655 
(1938). 
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TABLE II. Variation of sparking potential with spark 


56cm Vi kv Va kv V 


exp kv 

0.5 18.9 18.25 17.3 

1.0 32.8 32.2 31.6 

2.5 71.3 70.5 73.0 
5.0 132 132 138 
10 249 249 265 
15 361 363 386 


20 473 474 510 








Vi by »=10*%. Vy by Meek’s formula. Vexp: experimental values. 


DISCUSSION 


Thus as regards dependence of sparking 
potential in a uniform field both on spark length 
and on gas pressure there is little to choose 
between the values obtained by Meek’s formula 
and by the one based solely on the total number 
of ions produced in an electron avalanche during 
its passage across the spark gap. 

Considering the number of factors involved in 
spark formation it cannot be expected that any 
simple formula will accurately represent them 
all, even though the final solution is made to 
rest empirically upon the complicated experi- 
mental relation between a and X. Hence the 
discrepancies shown between the computed and 
the experimental results that have been discussed 
may not at first appear excessive. 

However, sparking potentials, owing to the 
rapid way in which a@ increases with X, are so 
extremely insensitive to very large changes in 
the exponential term e*®, which is the dominant 
quantity in each of the two formulas involved, 
that it is not justifiable to attach great signifi- 
cance to a hypothesis underlying any formula 
simply because this gives sparking potentials 
that are in fair agreement only with experimental 
values. This fact was recognized by Meek 
himself.’ 

The high degree of the insensitivity involved 
in the present problem is easily shown by using 
Meek’s formula backwards to find K for spark 
lengths other than 1 cm. Thus, assuming the 
correctness of the experimental results given in 
Table II for sparking potentials at different 
electrode distances, the corresponding values of 
X/p are first computed and from these a is 
obtained in each case from the experimental 
curve showing the relation between a and X/p. 
These two quantities are now substituted into 
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Meek’s formula to find the value of K for each 
gap. 

The values of K thus obtained for the sparking 
distances, 0.5, 1, 5, and 20 cm, for example, 
are 0.01, 0.09, 210°, and 2X10", respectively. 
Instead of being alike, as is postulated by 
Meek’s criterion, the two extreme numbers 
differ by a factor 210!’. In other words, when 
K is chosen so that the computed sparking 
potential agrees with its experimental value at 
5=1 cm, as was done by Meek, then the use of 
the same K for 6=20 cm gives a sparking 
potential only about 7 percent below that 
found by experiment, yet in order to bring this 
computed value into actual agreement with 
experiment it is necessary to increase K by a 
factor 2X10'*. Had Meek chosen the experi- 
mental value of the sparking potential for 6=5 
cm, for example, for determining the constant 
K in his formula he would have obtained for it a 
value 2X10, indicating an absurdly high radial 
field. It is thus by mere chance that K was 
found to be of the order of 1. 

From what has been said it appears that the 
results which have been presented lend no 
support whatever to the basic idea about radial 
fields on which Meek’s formula rests. 

When the relation n =e is tested in a similar 
way by finding the values of » which will make 
the computed potentials for different sparking 
distances agree exactly with those found by 
experiment, the disagreement among the values 
of nm is equally as large as that found above for 
the K’s in Meek’s formula. 

Since the exponential term e* is the all 
important quantity in both of the equations 
under consideration, it evidently does not 
represent properly the processes going on in the 
spark gap. This exponential term gives the total 
number of ions produced in an electron avalanche 
while it is traversing the whole gap on the 
supposition that the field remains unchanged. 
Actually, the space charges developed in the 
avalanche alter the applied field so markedly as 
to make calculations based on a uniform field 
altogether invalid. As has already been shown, 
these changes in the field are far from simple and 
it is not possible as yet to subject their effect 
to computation. 

It may well be that for voltages near the 
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sparking threshold, especially in case of long 
sparks, the initial avalanche does not continue 
to grow in intensity after it has built itself up 
to some limiting value, but owing to changes in 
the field is able subsequently to produce only a 
constant number of ions per centimeter of its 
path as it proceeds to the anode. After the 
electrons from the head of the avalanche have 
all entered the anode, the increase in the positive 
charge which results may then give rise to a 
positive streamer which in like manner makes its 
way to the cathode establishing a conducting path 
between the electrodes, as has been postulated. 

Once the positive streamer reaches the cathode 
it may not be necessary to look for processes 
other than the ones usually active at this 
electrode in continuous discharges through gases, 
for supplying the channel with electrons whose 
number will be limited only by resistance and 
other factors in the circuit. 

When the current in a spark has attained a 
very high value and a near equality between the 
electrons and positive ions in the channel has 
established itself, it is possible that the cross- 
sectional area of the channel, now at a high 
temperature, is prevented from expanding owing 





to the so-called pinch effect due to the magnetic 
field. 

Loeb and Meck'® have recently added a new 
process to the theory of the mechanism of a 
spark. They state that after a positive streamer 
has reached the cathode there a rush of 
electrons from the metal into the avalanche 
which gives rise to a high potential wave that 
passes up the channel to the anode and multiplies” 
the electrons there present by a large factor, 
thus making the channel highly conducting. 
I can find no evidence that the fields in waves of 
that character are sufficiently high to give 
electrons speeds adequate for the ionization of 
air at atmospherio pressure. 

Snoddy, Dietrich, and Beams"? measured the 
fields in such impulse waves started by a sudden 
flow of charge from a condenser charged to 
125 kv and found as the pressure of the gas was 
raised from 0.04 mm to 1.3 mm that the field 
changed from 2300 volts/cm to 420 volts/cm. 
Presumably at higher pressures the field would 
be still weaker. Such fields are far too small to 
produce ionization in air at atmospheric pressure. 

16 Reference 10, p. 40. 


7L. B. Snoddy, J. R. Dietrich, and J. W. Beams, 
Phys. Rev. 52, 739 (1937). 
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The Photographic Action of Electrons in the Range 
Between 40 and 212 Kilovolts 


R. F. BAKER, E. G. RAMBERG, AND J. HILLIER 


Research Laboratories, RCA Manufacturing Company, Camden, New Jersey 
(Received April 7, 1942) 


The Hurter and Driffield curves and the absolute sensitivities of Eastman medium lantern 
slide plates are determined for electrons accelerated through differences of potential from 40 to 
212 kilovolts. Within experimental error the H and D curves are identical in shape for the 
different voltages. The sensitivity reaches a peak near 100 kv and then drops off rapidly. This 
phenomenon may be utilized to increase the effective sensitivity of plates at high voltages by 
applying a thin, high density intensifying screen to the surface of the emulsion. 


INTRODUCTION 
HE advent of the electron microscope has 
increased the need for information regard- 
ing the action of electrons on photographic ma- 
terials. Thus, any translation of density varia- 
tions in an electron microscope picture into mass 
thickness variations of the object! must begin 


1 L. Marton and L. I. Schiff, J. App. Phys. 12, 759 (1941). 
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with a knowledge of the change in the photo- 
graphic density of the sensitive emulsion with the 
amount of charge per unit area which has reached 
the emulsion during exposure. Again, an acquaint- 
ance with the variation of the sensitivity and 
contrast properties of the emulsion used with the 
accelerating voltage employed is of importance in 
determining the most suitable operating voltage 
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and may lead to methods for enhancing the 
effective sensitivity of the recording materials. 

Earlier measurements of the electron sensitivity 
and of the density-exposure curves of photo- 
graphic materials have been carried out with 
electrons having velocities up to 100 ekv on the 
one hand,?~‘ and with beta-rays from radium B 
and C, with velocities in excess of 200 ekv, on the 
other.® In the first range a consistent increase of 
sensitivity with voltage was found, at least for 
low photographic densities. The converse was 
found to hold for the beta-rays. These results— 
and those of a number of other authors as well 
have been discussed in detail by Charlesby,® who 
derived a general formula for the slope variation 
of the density-vs.-log (exposure) curves (H and D 
curves) in terms of the grain size and the ab- 
sorption and excitation properties of the grains as 
well as the absorption of the emulsion as a whole. 

The reciprocity law, stating that the photo- 
graphic density depends only on the total charge 
per unit area falling on the emulsion (for a given 
voltage) and not on the length of exposure, is 
generally accepted for fast electrons (see, e.g., 
Becker and Kipphan*® and Bothe’). Deviations 
from this law, in the direction of lower densities 
for shorter exposures, have been noted by Borries 
and Knoll,‘ working with exposures of the order 
of 10-* second. It appears safe, however, to re- 
gard the reciprocity law as valid at least within 
the range of exposure times used in the electron 
microscope. 

The contrasting results of Schaffer,’ Becker and 
Kipphan,? and Borries and Knoll‘ for voltages 
less than 100 kilovolts, and of Aston and Ellis® 
for voltages above 200 kilovolts make the exist- 
ence of a maximum of photographic sensitivity in 
the region from 100 to 200 kilovolts a practical 
certainty. The presence of this maximum also 
follows from the commonly accepted laws which 
determine the energy loss of electrons in matter 
if the photographic density is assumed to be pro- 
portional to the energy per unit area absorbed by 





2 A. Becker and E. Kipphan, Ann. d. Physik [5], 10, 15 
(1931). 

3H. Schaffer, Arch. f. Elektrotechnik 26, 313 (1932). 

4B. v. Borries and M. Knoll, Physik. Zeits. 35, 279 
(1934). 

°C. D. Ellis and G. H. Aston, Proc. Roy. Soc. A119, 645 
(1928). 

6 A. Charlesby, Proc. Phys. Soc. 52, 657 (1940). 

7™W. Bothe, Zeits. f. Physik 8, 243 (1922). 
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the emulsion. These laws are: (1) the Lenard law 
of absorption that 


I= Ioe—**, (1) 


where Jo is the incident electron current and J 
that which remains at a depth x. a is a constant 
which, in the region up to about 50 kilovolts, ac- 
cording to Terrill’s measurements,’ is propor- 
tional to the inverse fourth power of the mean 
electron velocity and to the density of the ma- 
terial ; and (2) the Thomson-Whiddington law® 
governing the mean loss of kinetic energy eV of 
the electrons in passing through a layer of matter 
of thickness x. This law is stated as 


Ve— V*?=bx, (2) 


where the constant } is approximately propor- 
tional to the density of the material. 

The Thomson-Whiddington law by itself pre- 
dicts an almost linear increase of the energy loss 
in the emulsion (for constant incident current) up 
toa voltage V»= (bx), since up to this point prac- 
tically all of the electron energy is dissipated in 
the emulsion. Beyond this point an increasing 
proportion of the electrons (according to Eq. (1)) 
with an increasing amount of kinetic energy (Eq. 
(2)) pass through the emulsion into the backing. 
The maximum of energy loss within the emulsion 
occurs at the voltage for which the rate of in- 
crease in the total energy of the electrons leaving 
the emulsion is just equal to the rate of increase 
with voltage of the energy of the incident elec- 
trons. This maximum must exist if the absorption 
coefficient a decreases with increasing voltage in 
proportion to the reciprocal of the latter or more 
rapidly than this. This is true since, for Very large 
values of Vo, by Eqs. (1) and (2), the rate of 
change of the energy loss in the emulsion is given 


by 
d b = d(aVo) 
—(InVo= IV) = — Tox ~), (3) 
dVo 2Vo* = dV 





which becomes negative under the conditions 
stated. 


8H. M. Terrill, Phys. Rev. 24, 616 (1924). 

* For the constant 6, H. M. Terrill, Phys. Rev. 22, 101 
(1922), finds the value 4-10"p volt?/cm, p denoting the 
density, for the constant a, 12p/8*, 8 being the ratio of 
mean velocity of the electrons to the velocity of light 
(reference 8). 

1 For a different formulation of the energy-loss relation, 
see E. J. Williams, Proc. Roy. Soc. A130, 310 (1932). 
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Equations (1) and (2), however, are not ade- 
quate for deriving the variation in photographic 
density with voltage, for the following reasons: 
(1) the variation of a with voltage is not known 
beyond about 50 kilovolts; (2) the Lenard law 
fails at low electron velocities ;* (3) it is necessary 
to take account both of the velocity distribution" 
and of the angular scattering of the electrons in 
the emulsion, the latter determining the effective 
values of x in the Thomson-Whiddington law; 
and (4) there is some evidence that radiation 
emitted by the glass backing of photographic 
plates under the influence of electron bombard- 
ment may contribute to their blackening.‘ Even 
so, it is evident from the above considerations 
that the variation of photographic density with 
voltage must vary with emulsion thickness, the 
voltage corresponding to maximum photographic 
density for a given charge density increasing with 
the latter. 


MEASUREMENTS 


All of the measurements reported in this paper 
were made with the RCA high voltage electron 
microscope” with the object holder removed. 


"E. J. Williams, Proc. Roy. Soc. A125, 420 (1929). 
2V. K. Zworykin, J. Hillier, and A. W. Vance, J. App. 
Phys. 12, 738 (1941). 














Eastman lantern slide plates (medium) with the 
same emulsion number were used throughout, as 
these exhibit favorable characteristics with re- 
gard to grain size, sensitivity, and contrast prop- 
erties for electron microscope work. Development 
was standardized at 3 minutes at 19°C in 
Eastman D-72 developer diluted in a ratio 1 : 2 
with water. 

The measurements may conveniently be di- 
vided into two parts: the determination of the H 
and D curves for the different voltages, using 
relative values for the exposure ; and the determi- 
nation of the plate sensitivities in terms of the 
number of coulombs/cm?*, required to produce a 
given blackening. In order to obtain data for the 
H and D curves, the projector lens current of the 
electron microscope was adjusted so as to give a 
large magnification and the mask above the plate 
was adjusted to expose a single }” strip of the 
2”"X10” plate at a time. With the microscope 
voltage placed successively at 40, 78, 127, 170, 
and 212 kilovolts, strips were exposed for 5, 10, 
20, 40, and 80 seconds, leaving clear spaces be- 
tween adjacent steps. Nine such plates, with 25 
photometer strips apiece, were taken. The densi- 
ties of the strips, as compared with those of the 
intervening clear spaces, were measured with an 
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Fic. 2. Measurement of absolute plate sensitivity. 


improvised microphotometer consisting of an 
RCA 931 electron multiplier tube’ operated with 
a regulated power source and mounted on a plate 
with an aperture and inserted in place of a plate 
holder in the camera of a Leitz panphot micro- 
scope. The magnification used was between 60 
and 100. The output currents, ranging from 0.5 
to 400 microamperes, were measured with an 
ultrasensitive d.c. vacuum-tube microammeter." 

After conversion of the current readings into 
densities, the latter values plotted against the 
logarithm of the exposure yielded nine sets of five 
points for each voltage. Shifting the sets of points 
along the exposure axis so as to cause them to 
fall, as nearly as possible, on the same curve in 
order to account for differences in the beam in- 
tensity used for the different plates, it was pos- 
sible to obtain an average H and D curve charac- 
teristic for each voltage. The curves so obtained 
were found to be identical within experimental 

13]. A. Rajchman and R. L. Snyder, Electronics (De- 


cember 1940), p. 20. 
14 A.W. Vance, Rev. Sci. Inst. 7, 489 (1936). 
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error. Bringing the five curves for different volt- 
ages, again, as nearly as possible to coincidence 
by shifts along the exposure axis resulted in ob- 
taining the mean H and D curve shown, together 
with all the individual measured points, in Fig. 1. 

The measurement of the absolute plate sensi- 
tivities requires a measurement of the current 
producing the blackening on the plate. For this 
purpose a Faraday cage in the form of a nickel 
cylinder 7; inches in diameter and 13 inches in 
length provided with an aperture 0.150 inch in 
diameter was placed above the shutter (Fig. 2), 
and centered so as to receive all the current 
passed by the limiting aperture A above the pro- 
jector lens when the lens current was reduced 
to a very small value. Quartz insulation was pro- 
vided for the lead from the cage, which was car- 
ried out through a seal in a Pyrex window of the 
viewing chamber and connected through a shielded 
cable to the ultrasensitive d.c. vacuum microam- 
meter. The currents measured were of the order 
of 5-10-'° ampere; a bias of +90 volts did not 
affect the reading. 

After the meter was read, the Faraday cage 
was swung out of the way with the aid of a hook 
mounted on a sylphon bellows, the projector lens 
current was increased until there was a luminous 
disk approximately 4 cm in diameter on the 
fluorescent screen, and a 10-second exposure was 
made. After development of the negative, carried 
out as before, each ‘disk’? was marked with 
equally separated rings along two diameters, as 
shown in Fig. 3a. Density measurements were 
made in each ring and the results for rings lying 
on a circle about the center averaged, yielding a 
curve for the variation of the density outward, 
D(p) (Fig. 3b). This was found necessary, as even 
a small amount of distortion of the projector lens 
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Fic. 3. Density distribution on plate: (a) Marking of plate; 
(b) example of density distribution. 
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results in a considerable dropping off of intensity 
in the marginal regions of the exposed area. 
Since it had been found that the H and D 
curves for different voltages were the same except 
for a shift along the log (exposure) axis, the true 
exposures in coulombs/cm* at the voltage V re- 
quired to produce the densities D are given by 


Ey(D)=E(D)\Cy, (4) 


where E(D) is given by the curve in Fig. 1. The 
constant Cy, on the other hand, is a direct meas- 
ure of the relative sensitivity of the plate at 
different voltages. Now the total amount of 
charge reaching the plate during exposure is 


Pi 
an f Ev(p)pdp=1-t, 
0 


where p; is the radius of the blackened circle, 7 
the measured current, and ¢ the exposure time. 
Since 


Ey(p) = E(D[p }) ‘Cy= E(p) ‘Cy, 


where E(p) can be obtained from D(p) with the 
aid of the curve in Fig. 1, Cy can be determined 
by a graphical integration : 


(5) 


Figure 4 shows the results of four measurements 
of Cy for each of the voltages 40, 79, 127, 170, 
and 212 kilovolts and of two measurements for 








the intermediate voltages 67, 96, and 105 kilo- 
volts, together with their probable errors, the 
latter being indicated by the length of the line 
segments. 


DISCUSSION AND APPLICATION 


Figure 5 shows the H and D curves for Eastman 
lantern slide plates (medium) in absolute units, 
as determined from the above data, and compares 
them with the results of Becker and Kipphan? on 
Agfa Extrarapid plates and of Borries and Knoll* 
on Agfa Isochromfilm, the exposure times in the 
latter case being of the order of 10~-* second. As 
was to be expected, the lantern slide plates are 
somewhat slower than the Extrarapid plate. 
Apart from this the principal difference between 
the present results and those of Becker and 
Kipphan is the absence of intersections between 
the curves in the present measurements. 

The fact that for voltages in excess of 100 
kilovolts the plate sensitivity rises rapidly with 
decreasing voltage may be utilized to enhance the 
effective sensitivity of the plate at high voltages 
by means of an “‘intensifying screen.”’ If this is in 
the form of a thin film of high density (e.g., gold 
or silver), the velocity reduction as well as the 
angular scattering within the film will enhance 
the action of the electrons on the emulsion with- 
out appreciably affecting the sharpness of the 
image. 

Figure 6 shows the effect of an evaporated 
silver film with a mass thickness of 3.4 mg/cm? 





mi EASTMAN LANTERN SLIDE, MEDIUM 
cou © w-oqee--. LANTERN SLIDE, COATED WITH 
70 CELLOPHANE, 3.3mg/cm* 
2 .--{g----- LANTERN SLIDE, COATED WITH 
7? CELLULOSE SCOTCH TAPE, 9.6mg/cm* 
co single a \ 
4 Measurement ‘ 
2.5). . 
3 \ ‘. 
2 \ ‘, a —— 
a . 
: \ 
20 5 _ er fl * < — . 
: 4 \ : seal Fic. 4. Variation of 
F plate sensitivity with ac- 
5 St celerating voltage 
use 3 yj tii Psciaied « c . 
° 
7) 
° 
5 
. 
os. ¥ 
— 
a 
“. 
bad . 
° a o 1 ai 1 i l 1 l 1 i 
20 40 60 60 100 120 140 160 180 200 220 240 ts KV 
Accelerating voltege corresponding to electron velocity 
454 JOURNAL OF APPLIED PHYSICS 





CS 





D (density ) 


18 


—  §ASTMAN LANTERN SLIDE, MEDIUM 


anc ~~ <deiimmbia AGFA EXTRARAPID (BECKER AND KIPPHAN [2)) 


‘ 
08 ‘ 
s 


06) i = s 





—-—— AGFA ISOCHROMFILM (BORRIES AND KNOLL [4] ) 








12.0 12.2 12.4 12.6 128 41.0 12 11.4 


‘an re zs 
o2— ae - peo 
— - ee 
Oe } 
I - i Se Se l i l 1 { L. | 4 ] I ! 
3 11.6 1.8 10.0 10.2 10.4 10.6 10.8 


fog, Ey (coulombs/cm* ) 


Fic. 5. H and D curves for different accelerating voltages and different negative materials. 


placed directly on the emulsion at a beam voltage 
of 127 kv. The film covered the area with the 
lighter background in the photograph. It is seen 
that the sharpness of the image (zinc oxide 
smoke) does not suffer due to the presence of the 
film where the latter adheres closely, even though 
the reproduction represents a sixfold optical 
enlargement of the original negative. At the edge 
of the film there results a brightening of the back- 
ground and a darker appearance of the film itself 
from the lateral scattering at the (slightly raised) 
film. 

In order to obtain a quantitative estimate of 
the effect of such ‘intensifying screens’ on the 
sensitivity of the plate at different voltages, 
Cellophane foil and cellulose Scotch tape were 
used to coat a portion of a lantern slide plate 
exposed to the unobstructed electron beam in the 
electron microscope, the voltage being varied 
from 40 to 212 kilovolts. The use of cellulose in 
place of the more suitable heavy metal screens 
for intensity measurements is justified since the 
degree of scattering and absorption of high speed 
electrons depends primarily only on the mass 
thickness of the scattering and absorbing layer. 

After standard development of the negative 
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the densities of the portion exposed under the 
screen and of that exposed directly were meas- 
ured. If the same shape is assumed for the H and 
D curve for the plate with the intensity screen as 
without, the difference in log E corresponding to 
the two densities, according to the general H and 
D curve in Fig. 1, yields the ratio of the exposures 
required to produce the same density on the 





Fic. 6. Electron micrograph of zinc oxide smoke taken 
on plate partly covered by silver film (3.4 mg/cm?). Mag- 
nification optical enlargement 6X. 
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coated and the uncoated plate. The validity of 
assuming the same shape of the H and D curve 
for the coated as for the uncoated plates was 
established by making separate measurements of 
the curve for a plate coated with cellulose Scotch 
tape at 127 and 170 kilovolts, i.e., at voltages for 
which the screen definitely reduces and increases, 
respectively, the sensitivity of the plate. As 
shown in Fig. 7, the shape of the H and D curve 
for the coated material does not deviate from the 
general curve for the uncoated plate by more than 
experimental error. Multiplying the sensitivity 
curve (Fig. 4) for the uncoated plate by the ex- 
posure ratios corresponding to the differences in 
density of the coated and the uncoated plate 
leads to the sensitivity curves for the coated 
plates, as shown by the dotted lines in Fig. 4. It 
is seen that 

1. At low voltages the sensitivity is reduced due to the 

absorption and energy losses of the electrons in the 
screen. 

. The maximum of sensitivity is shifted toward higher 

voltages in proportion as the screen is made thicker. 

3. For thin screens the sensitivity maximum is higher than 
for the uncoated plate, due to the lateral scattering 
within the screen, causing more oblique incidence of the 
electrons on the emulsion. 


Nm 


4. For thicker screens the maximum is lower and broader, 
due to appreciable absorption in the screen and to the 
large velocity spread of the transmitted electrons. 
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An apparently simpler solution for utilizing the 
intensity maximum near 100 kilovoltsin recording 
high voltage electrons consists in introducing 
decelerating electrodes which bring the kinetic 
energy of the electrons down to about 100 kv 
before they reach the plate—just as, in cathode- 
ray oscillographs, electrons are at times given a 
final acceleration in order to increase the light 
output of the fluorescent screen. More often than 
not, however, the introduction of the necessary 
additional high voltage electrodes will constitute a 
serious drawback for construction and operation. 

Again, an even somewhat greater increase of 
the plate sensitivity at high voltages than with 
the nassive intensifying screen could be attained 
by increasing the emulsion thickness by a corre- 
sponding number of mg/cm?. However, the added 
thickness would, due to the smaller density of the 
emulsion, have to be larger than for a heavy 
metal screen, the spreading of the electrons in the 
emulsion would be greater, and correspondingly 
the resolving power of the plate would be reduced 
by a larger amount. 

We wish to express our appreciation of the 
many valuable suggestions given by Dr. V. K. 
Zworykin, under whose direction this work was 
carried out. 
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The Torsion of a Prism with Cross Section the Inverse of an Ellipse 


THOMAS JAMES HIGGINS 
Tulane University of Louisiana, New Orleans, Louisiana 


(Received May 1, 1942) 


Saint-Venant’s torsion problem is solved for the linear prism with right cross section bounded 
by the curve, an inverse of an ellipse with respect to its center. By a method developed recently 
by R. M. Morris, explicit formulas are obtained for the torsion, conjugate torsion, and complex 
torsion functions; for the twisting moment; for the stress function and the components of 
shearing stress. This work is of interest as complementing recent work by A. C. Stevenson and 
by D. L. Holl and D. H. Rock; further, as an addition to the very few exact solutions that have 
been effected for cross sections bounded by single, closed, continuous, partially concave curves. 


A C. STEVENSON! and D. L. Holl and D. H. Rock? published, independently and nearly 
¢ simultaneously, solutions of Saint-Venant’s torsion problem for the linear prism with right 
cross section bounded by an elliptic limacon. This curve is the inverse of an ellipse with respect to a 
focus. As complementing the work of these three it is of interest to possess the solution for a kindred 
curve, the inverse of an ellipse with respect to its center—i.e., the curve v=0 in the net defined by 


z=csec(wt+ik), 2=x+1y, w=u-+1, (1) 


c and k constants defining the semi-major and semi-minor axes of the ellipse of inversion, respectively 
directed along the x and the y axes and of length c~ cosh k and c~ sinh k. Additional interest attaches 
to this contour in that a recent exhaustive survey of the known analytical solutions of Saint-Venant’s 
torsion problem (by the writer; to appear elsewhere shortly) revealed that exact solutions have been 
effected for but very few cross sections bounded by curves of this type; that is, by single, closed, 
continuous curves possessing concavity. 

We utilize the powerful method of treating two-dimensional physical problems recently developed 
by R. M. Morris.*~> With respect to the torsion problem for singly-connected cross sections, Morris 
has proved that if the boundary curve is v=0 in the net defined by 


@ 
z= > ane’™, s=x+iy, w=u+in, (2) 
m=0 
and if v= is the sole singularity within the cross section, the complex torsion function is 


Q=ibo/2+i ¥ dye, (3) 


m=1 


and the twisting moment is the real part of 


ea oa 
M= (rum/2)[cobo+ > a (CmOm + C—mDm) —2 z. MDnbm |, (4) 
m=1 m=1 
where 
«© ™ oa 
bm= Do GAm+s 5 bm= LS Ams, (5) 
s=0 s=0 
x ea 
Cm= DL (M+S)Admye3 Cm= LDL Sdddimys, (6) 
s=0 s=0 


the bar indicating the complex conjugate. 





1 A. C. Stevenson, Proc. Math. Soc. London [2] 45, 126 (1939). 

2D. L. Holl and D. H. Rock, Zeits. f. angew. Math. u. Mechanik 19, 141 (1939). 
> R. M. Morris, Math. Ann. 116, 374 (1939). 

4R. M. Morris, Math. Ann. 117, 31 (1939). 

®R. M. Morris, Proc. Math. Soc. London [2] 46, 81 (1940). 
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It is convenient to note here the following identities: 


sec (w+ik) =2 x (—1)0—D/2errkgirw (r=1,3,5---), (7) 

csch k= 2 De™ (y=1,3,5---), (8) 

esch (2k) coth (2k) =2 rr -2rk (r=1,3,5---), (9) 
esch? 2k=2 E ne—2nk (n=2,4,6---), (10) 

1—coth k= —2 5 2m (n=2,4,6---). (11) 


n=2 


In light of these identities we understand hereafter that r designates only an odd integer, and that n 
designates only zero or an even integer. 
Use of (7) enables us to write (1) in the form of (2) 


® 
2=2c > (—1)¢-» 2e—rkeirw . (12) 
r=1 
whence comparison of (12) with (2) reveals that 
dm=0, m even or zero; dm=a,=2¢c(—1)°—/*e-"*, m odd. (13) 
By (13) a,4,4, vanishes except if s is odd and m is even. Accordingly, by (5) and (13) 
™ 4) 2) 
bm = Om =bn= FY Adnyr=4e7(—1)"/2e-"* DO et; (14) 
r=1 r=1 


whence by (8) 
b,, = 2c? csch (2k)(—1)"/*e-™. (15) 
Substituting (15) in (3) and utilizing (7) yields 
Q=ic? csch (2k)[1—e'”*+™ sec (w+ik) | 
=c* csch (2k) tan (w+7k). a 


The real and imaginary parts of Q are, respectively, the torsion function ¢ and 7 times the conjugate 
torsion function y. Utilizing the relationships ¢=(Q+9)/2 and y=(Q—Q)/21, simple trigonometric 
analysis yields 


Q=o+i~=c? csch (2k)[sin 2u+i sinh 2(v+k) ][cos 2u+cosh 2(v+k) ]-. (17) 
The stress function V is defined by 
VW =y—(x?+ y*)/2. (18) 
As by (1) 
x?-+-y? = 22 =c* sec (w+ik) sec (W—ik) = 2c*[cos 2u+cosh 2(v+k) |-', (19) 
we have from (17) and (18) 
v= —c*[1—csch (2k) sinh 2(v+k) |[ cos 2u-+-cosh 2(v+h) }-'. (20) 
The shearing stresses are given by 
X,=TpdWV/dy; Y,=—rpdV/dx. (21) 
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we obtain by carrying out the simple but lengthy analysis indicated 

















Tesch (2k) cosh (v+k) sinh (v+k) 7 

X,= —2rwe sin u —- —|, (24) 
Lcos 2u—cosh 2(v+k) cos 2u+cosh 2(v+) J 
resch (2k) sinh (v+k) cosh (v+k) , 

Y,= —2rpc cos u _ , (25) 
|cos 2u—cosh 2(v+k) cos 2u+cosh 2(v+k). 


Calculation of the twisting’ moment requires evaluation of the ¢ coefficients. Noting the remark 
preceding (14), we have from (6) and (13) 


Coim=Con= D> 1 npe=407(—1)*/2%e-** > re 2k, (26) 
~~ r=] 
whence by (9) 
C_»=2c? csch (2k) coth (2k)(—1)"/*e-™. (27) 


The a coefficients being real, we have by (5) and (6) 


oo 


Cm=Cao= D> (nt) ddayr=nbnt+c_n. (28) 

r=1 

Accordingly, in view of (28) we have for (4) 
M=(rux/2)[cobo+ X (2¢-nba—nbn2) (29) 


Evaluating (29) by use of (15) and (28) yields 


M = (ryum/2)[4c* csch? (2k) coth 2k+8c! csch? (2k) coth 2k + 


9 


«2k 4c4 csch? 2k } ne~?"*], (30) 


whence by (10) and (11) 
M =rurc*(2 csch? 2k+csch‘ 2k). (31) 


As the semi-minor and semi-major axes of the cross section are, respectively, a=c sech k and 
b=c csch k, (31) can be transformed easily into 


M = tur(a*+b'+6a7b*) /16. (32) 


Corroboratively, if a=b, the original ellipse is a circle, the inverse with respect to its center is a 
circle, and (32) yields the known expression for the twisting moment of a torsed circular cylinder: 
M =ryra'‘/2. 

From (24) and (25) the shearing stresses on the contour at the ends of the semi-major and semi- 
minor axes are 


X ,= —7tp(a?+5b?)/2b; YV,=rp(a?+5?)/2a. (33) 


Corroboratively, if a=b, (33) yields the known expression rya for the shearing stress on the surface 
of a torsed circular cylinder. 
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Human Trajectories 


RuFUs OLDENBURGER 
Illinois Institute of Technology, Chicago, Illinois 
(Received May 21, 1942) 


It is proved that a human being M falling from a ledge A through the air to the ground 
below will attain the greatest range if M falls forward while standing upright. If M was found 
on the ground at a position B beyond this range of fall, M furnished energy himself to attain B. 
An analysis of trajectories from A to the ground due to an effort on the part of M is made. 
If the kinetic energy which ./ must produce to reach B is much greater than that which WM 
can furnish after loss of balance, M made the trip from A to B by a premeditated jump. Such 
analyses are important for determining suicides. 


INTRODUCTION in the air the center of mass C moves approxi- 


STUDY of the possible motions of a human 

being M from an initial position A to a 
final position B will be made here by the use of 
the fundamental laws of motion only. 


mately along a parabola, we may restrict our 
study to plane motions. Because on contact with 
the ground the center of mass C no longer con- 
tinues to move in a parabola, we shall analyze 





We let M stand upright on a horizontal rod R. y 
We let / designate the distance from the center 
C of mass of M to R and let r represent the radius 
of gyration of M with respect to R. For a par- 
ticular body M, the numbers / and r can be 
obtained by the usual balancing and pendulum 
methods. 

We suppose that air resistance can be neg- 
lected,* as is true, for example, if B is no more 
than 100 feet below A. Since once the body M is 








y 


? 


a 


FIG. 
our motions as if C continued to follow the 
parabola traversed by C before contact with the 
ground. In this case C would strike the ground at 
a point D. We let b designate the distance A is 
above D. We set up a Cartesian coordinate 
(x, y)-system with origin at A, the x-axis hori- 
zontal, and D to the right of A. We let 
D*=(a, —b) be the coordinates of the given 
point which corresponds to the actual motion 
taken by the body M from A to the horizontal 
line y= —b, which motion is to be determined 
from A and D*, as well as other properties of the 
motion. (See Fig. 1.) 


Be 


W 





——ge XX 


Cc 
y--£ 





ANALYSIS OF THEORY 
1. Falls 





Mc-4) 2*(a,-4) 


Fic. 1. 


We shall consider a fall from A to the line 


*Computation shows that a 10 m.p.h. wind would . . 
; y=-—b. We let (c, —b) denote the coordinate of 


deviate M by less than an inch in a drop of 30 feet. 
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Fic. 3. 


the point D in this case. We shall determine the 
maximum value of c for all possible falls from 
A to the line y= —b. 

We suppose that / stands erect at A and falls 
forward. Until M departs from the point A the 
distance AC=h remains practically unchanged. 
We let 6 be the angle measured from the positive 
half of the y-axis to the line AC at any time before 
departure of M from A (Fig. 2). We let m and g 
designate the mass and acceleration of gravity 
of M in the usual units. The torque L on M 
with respect to A is mgh sin 6. We let J denote 
the moment of inertia of M with respect to A, 
whence J=mr*. From the fundamental relation 
L=I]60"', where 6” is the angular acceleration of 
C, we have 
- sin @. 
r2 
Integrating, using the initial condition that 6’=0 
when @=0, we obtain 


2(gh)} sin 6/2 
Ol tive, (1) 
r 
The centrifugal force F, acting on M is 
F.=m(6’)*h. 


When this force equals the component ga. of 
gravity along the line AC, the body M flies out 
into space, the constraint being zero. Equating 
F. and gac we obtain 


2mgh? 
mg cos 6@=— 





: (1—cos @), 


hd 
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whence 
2 


2+r? 


where r=7r/h. The solution @ of (2) will be called 
the ‘‘angle of departure.”” Experiment shows that 
for a body M, the number 7? is about 3, whence 
6=55°. If we let 0) designate the solution of (2) 
(see Figs. 3 and 4), by (1) and |v| =/6’ the linear 
velocity v of the point C at the instant fo of de- 
parture has the magnitude 


cos 6= 


(2) 


2(gh)' 60 
|= sin —, (3) 
2 


9) 


}¢U | 





y 


and is directed downward along the tangent to 
x?+y?=h? at. the point Co=(hsin 0, h cos 4) 
occupied by C at the instant to. The components 
v, and v, of v in the x and y directions have the 
magnitudes |v| cos 4, |v| sin 4, respectively. 
Upon departure the point C travels along the 
well-known parabola 


y—k=(x—J) tan c———___—— (4) 


where (j, k) =Co, a= —O6 (see Fig. 2). Setting 





Fic. 4. Falling rod which was released from the vertical 
position. Exposure frequency—15 per second. In 4th 
position, the rod makes an angle of 60° with the vertical, 
and has already left the ledge. 
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Q 
Fic. 5. 
y=-—b in (4), and solving for x we obtain the 


abscissa c of D. For @)=55° this abscissa c is 
given by 


c= .0025[5.74(56h— |v|2) 
+ |v| (33|v|*+448(7b+4h))!]. (5) 


In the above paragraph we supposed that M 
fell forward while erect. It remains to show that 
the range c thus obtained is a maximum for all 
falls from A to the line y= —b. We shall accom- 
plish this proof by showing that as the center of 
gravity of M is lowered, as when M crouches, 
the range on y= —b is shortened, i.e., the abscissa 
of the point where C meets y= —) is decreased. 
We shall actually carry through the converse 
argument by treating the trajectories as M 
straightens up. 

Experiment shows that when AC makes an 
angle of about 10° with the vertical, a body M 
loses balance and falls over. We shall term this 
angle the angle of loss of balance. It follows that 
M cannot lean over very far from the vertical 
without starting to fall over. In view of the con- 
tinuity of the functions involved, mathematical 
computations show that the substitution of a 
position of M where AC makes an angle 
6(0=@=10°) with the vertical for the position 
where AC is vertical will affect the range on the 
line y= —b by a relatively small amount. It will 
therefore be no restriction on the generality of 
the method to compare trajectories where AC is 
initially vertical. 

The trajectory of M is influenced most by the 
trunk and head of the body, and only in a small 
way by the legs. We shall, for the moment, 
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neglect the legs, and compare trajectories which 
arise when the trunk T of M is raised in such a 
manner that one vertical chord LN of M is kept 
vertical. We understand the trunk T to refer to 
the part of M above the legs. (See Fig. 5.) The 
center of gravity C is near the y axis, for which 
reason we shall suppose that C is on the y axis. 
We shall suppose that JT is raised through a 
distance e. It follows that the distance h from A 
to the center of gravity C increases from h to 
h=h+e. We shall determine how the radius of 
gyration r of T with respect to A varies under 
this elevation of T. The point (x, y) in T goes 
into (x, y+e) under this change. Thus the mo- 
ment J, of inertia of T with respect to the 
x axis goes into 


I,’ =1,+2em,+em, 


where m, is the first moment of T with respect 
to the x axis before T is raised, and m is the mass 
of T. Under the elevation of T the moment of 
inertia of JT with respect to the y axis remains 
invariant. Thus the moment of inertia J’ of T 
with respect to A after elevation is related to the 
corresponding initial moment of inertia J by the 
formula 


I’=I1,/+I,/ =I+2en,+em. (6) 
We let # designate the radius of gyration of T 
after elevation. By (6) we have 
where 
k=2eh+e. 


We also have h?=h?+k, whence 
= “ =——-, (7) 


where ph?=k. Since r>1 for M, we have 7<r, 
from which (2) implies that the angle 6 of de- 
parture for the elevated trunk is smaller than 
the initial departure angle 4). We can write 


(3) as 
2gh \3 
|v | -(—) : 
2+7? 


It follows that the magnitude of the new velocity 
vector 3 of T is greater than |v| for the initial 
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departure velocity v ot T. The abscissa 7 of the 
point C of departure of T for e=0 is hsin 6, 
whence , 

_ h(o?+4e)! 


jp eee, (8) 
2+0o 





For small values of e, we have p<1, whence we 
can write ¢= 7", given in (7), as 
o+p(1—o)(1—p+p*—---). 


Writing Ao=é¢—«¢, Ah=e, expressing p in terms 
of Ah, and observing that p—0 as e—0, we have 


Now (8) implies 


dj o*(o+6)+8 


— = 9) 
dh (o?+4oe)}(2+6)? 





Certainly 
dj/dh>0, (10) 


whence j increases with an increase in h. Thus 
the abscissa 7 of the point (7, k) of departure of 


C, where AC=ii+e=h is greater than j.* It 
follows that we have the situation illustrated in 
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Fig. 6. We have proved the validity of the set 
of inequalities |5|> |v], 7>7, h>h, 60<6. If 
we replace the trajectory S determined by @ 
at P by a trajectory Q determined by 6 acting 
at P, the range on y= —b is decreased. Thus, 
if the range of Q on y= —b exceeds the range of 
the trajectory H, determined by v at P, on 
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y=-—b, this is certainly true if Q is replaced 
by S. Since we have reduced the comparison of 
S and H to that of trajectories Q and H emanat- 
ing from the same point P, we shall compare 
corresponding trajectories Q’ and H’ passing 
through O with the velocities 6 and v, respec- 
tively. Our equation of motion is (4) with 
(j, k) =(0, 0). The slope dy/dx of the trajectory 
(4) through A is given by 


32x 





y’ = tan a——————_. 
\v|* cos? a@ 


We assign to x a positive value K. We are in- 
terested in negative values of a. Now the absolute 
value | y’| of y’ is 


—tan a+32K/|v\* cos® a. 


As |v| is increased, | y’| is decreased, whereas as 
'a| is decreased, |y’| is also decreased. We let 
Qx and Hx designate the points on Q’ and IH’, 
respectively, obtained by setting x=K (Fig. 7). 
The slope of Q’ at Qx is numerically less than 
that of H’ at Hx. Since v makes a larger angle 
with the x axis than 3, the point Qx is above Hx 
for small values of K. On the other hand, for Q’ 
to cross H’ at some point Qx=Hx, the slope of 
Q’ must be greater numerically than that of /7’ 
at this point. It follows that the trajectory Q’ 
lies above H’ for each K. Thus the range of Q’ 
on each line y= —6 is greater than that of #7’. 
It follows that the range of S on y=—6 is 
greater than that of Q on this same line, whence 
elevating T increases the range on y= —b. 

In the above treatment of JT, we used the 
property that the motion of M is dominated 
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by 7. We shall now show that the effect of the 
legs on the trajectory leaves invariant the prop- 


erty that the range on y= —d is increased as T 


is raised. When M is raised so that T is elevated 
through a distance e while the feet of / are kept 
in contact with A, this change of the position 
of M may be thought of as being accomplished in 
two steps, namely, the elevation of MV as a rigid 
body though a distance ¢, and the subsequent 
dropping of the bottoms of the feet through a 
distance « as illustrated in Fig. 8. We imagine 
that a rigid weightless rod R is affixed to M and 


joins A in II along the y axis with the body M. 
We let II be the trajectory traversed by C as M 
affixed to R falls forward. We let II; and I be 
the corresponding trajectories for I and III, re- 
spectively. In view of the development above 
the range of Ilj; on y= —d is greater than that 
of Il. We shall show that as we go from II to III 
we do not decrease the range of Iq; on y= —)d in 
such a way that this range for IIy1 is less than 
that of II;. We let h and 7 be the distance from 
A in III to the center of gravity of M, and the 
radius of gyration of M with respect to A, re- 
spectively. We let 4 and r be the corresponding 
numbers for I. We let (x, y) be the coordinates 
of a point P in the legs of M in the position I, 
and (#, y+e—’) the coordinates of P in the posi- 
tion III. We write 


. 


mé= | e'dm, 
l 


where m is the mass of M, and the integral is 
taken over the legs / of WZ. Since the trunk T of 
M is kept fixed in the transition from II to III, 
the first moment My, of T with respect to the 
x axis AX is invariant. The first moment M’,x 
of M in III with respect to AX is related to 





Fic. 9. These figures (20 exposures per second) show that the range of the center of gravity (knob) is increased as a 
weight is moved out on a rod. 
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this moment Myx for the position I by the 
formula 


M' 4x = Max—mi+me. 
Hence we have 
h=h+(e—@) (11) 
approximately, since the abscissa of the center of 
gravity C of M in III is very small relative to 
e—é. The moment of inertia I’4y of M in III 
with respect to AX is related to this moment of 
M in I by the formula 
TI’ ax =Iax+[2(eh—p)+(C+0*) Jm, (12) 


where 


um= f (y+ededm, v'm= { (¢)*dm. 
1 l 


Thus approximately (neglecting infinitesimals of 


higher order than 1) by (11) and (12) 


r?+2eh—2y 
#2 = —___ (13) 


For é small relative to e, we have 
(1—é/e)r>h. 


For a body M this last inequality is satisfied, 
whence 7<+r. Thus the angle 69’ of departure for 
III is less than the corresponding angle 4 for I. 
Now the velocity of departure v’ for III has a 
magnitude greater than that for I, since +<r, 
and h>h. We let 7; be the number for II corre- 
sponding to 7 for I. The ratio 7;? is approximately 


r?+2eh 
h?+2eh 


9 
Ty = 


(14) 


When h>y-+e for each point (x, y) in 1, we have 
éh>y. Now 1: is obtained from 7 in (13) by 
adding 2éh, and 24 to the denominator and 
numerator, respectively. Thus 


T>T>71. (15) 


In view of (15), the derivative do/dh{o = 7? | for 
changes of the type from I to III is numerically 
less than do/dh for changes of the type from I 
to II. It follows that the derivative of j as given 
in (8) with respect to h is at least as great for 
changes from I to III as for changes from I to IT. 
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Thus dj/dh satisfies (10). It follows from the 
theory developed above for transitions from I to 
II that the range of M for III on y= —b is at 
least as great as the corresponding range for I. 
We have proved that the effect of the legs is not 
enough to prevent the range from increasing as 
the trunk of M is elevated vertically. (See Figs. 
9 and 10.) 





Fic. 10. Comparing Figs. 9 and 10 shows that the effect 
of an additional weight on the rod near the ledge affects 
the trajectory of the center of gravity only negligibly. 
Exposure frequency for Fig. 10—15 per second. 


We have omitted above the treatment of the 
case where the axis of the trunk T of M is in- 
clined from the vertical. We consider a trunk T 
which is symmetrical with respect to two mu- 
tually perpendicular lines CJ and CK through 
the center of gravity C of T. We shall compare 
the trajectories of T for an angle of 0° between 
CJ and the positive half of the y axis, where C 
is on the y axis, and the trajectory obtained by 
inclining CJ through an angle a about C as 
illustrated in Fig. 11. A point P with polar 
coordinates (p, ¥) in Fig. 11, I, relative to CK as 
axis, goes into a point with coordinates (p, ~—a) 
in Fig. 11, Il. The center of gravity of T remains 
invariant under transition from I to II. This is 
also true of the moment of inertia with respect 
to A, since the moment J, of inertia of T in the 
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position I with respect to A is given by 
La= | (i? +92-+2plt sin y)dm, 
T 


(dm is the usual element of mass) whereas this 
moment J,’ for II is 


La'=1s+2h f ofsin (y—a) —sin yp jdm, 
> 


whence using the symmetry properties in splitting 
the integral over 7 into integrals over the four 
regions obtained from T by using CK and CJ as 
cutting lines, we obtain 74’ =J,. Since the radius 
of gyration r of T with respect to A, and h are 





invariant, the trajectories obtained by letting MW 
fall forward, with 7 in the initial positions I and 
II are identical. (See Fig. 12.) As in the treat- 
ment of the elevation of 7, introducing the effect 
of the legs will cause only a small deviation in 
the range on y= —b. Also a small deviation of C 
from AY in the initial position of M, and a 
reasonable departure of JT from symmetry with 
respect to CJ and CK will have only small 
effect on the range y= —b as some simple trials 
of deviations show. 

If the body M falls forward while M is not 
resting on his feet, the center of mass of M is 
much closer to A than when WV is standing erect 
on his feet, whence, by simple mathematical 
considerations, the range of WM on y=—O6 will 
be considerably less than when .V/ falls forward 
from the erect standing position. It follows from 
all of these considerations that the most efficient 
fall of MW from A is attained when the initial 
state of M/ is the erect standing position. (See 
Fig. 13.) 

That a body M should fall in such a manner 
that the maximum range C in (5) is attained is 
exceedingly unlikely. If in a given physical case 


the range exceeds the C of (5) we may conclude 
that the body M furnished energy himself to 
attain the given range. 





Fic. 12. The trajectories of the centerZof mass for these figures are identical. Exposure frequency—20 per second. 
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Fic. 13. Experimental set-up for photographing falling 
rod. The author is shown on the left and Mr. Howard 
T. Betz of the Armour Research Foundation on the right. 
The above photographs were taken by Mr. Betz. 


From the well-known formula s= 16/7, we find 
that the angle the “axis” of MW makes with the 
positive horizontal direction, measured clock- 
wise, is by (1) given as 


sin 09/2 
69 +———[gh(b+h cos 6») |} (16) 
2r 


at the instant C strikes y= —b, provided the 
angular velocity of M has been constant during 
flight, a state which may be considered normal. 
From (16) one can determine whether or not a 
most efficient fall of M4 would land M on his 
head, or on his feet, or flat on the ground. 


2. Jumps 


In making a jump from the point A the 
center of gravity C of M must be within a dis- 
tance h of A at the instant ¢ of departure, 
whence C is within or on the circle [ with equa- 
tion x*+y?=/h? at this instant. We shall now 
determine the point which must be occupied by 
C if M is to leave A with the minimum speed, 
and thus minimum kinetic energy, and still 
reach D*=(a, —b), and we shall derive this 
minimum velocity and energy. 

The argument here was worked out by D. B. 
Dekker.f Solving (4) for |v| after setting (x, y) 


+ Assistant in the Department of Mathematics, Illinois 
Institute of Technology. 
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=(a, —b) we obtain 
16(a—j)* sec a 


[(a—j) tan a+(k+b)] 


9) 
)¢ 


We shall minimize |v| by varying a. Now |v| is 
a minimum when |v|? is a minimum. Setting 
d\v\?/da equal to zero, and solving we obtain 
—(b+k)+[(b+k)?+(a—j)?}} 
tan ¢ 2 —————————————. (97) 
(a—J) 

We have used the + sign in (17), because M 
cannot jump downward, and the — sign would 
by j<a yield a negative tan a and negative a. 
For a@ as given in (17) we have 


|v|*=32(LO+e)*+(a—j*P—(+k)}. (18) 


At the point (j,k), the speed |v| of C as given 
in (18) is the least speed with which C can depart 
from (j,k) and reach D*. At a particular point 
(j,k), the speed |v| of (18) takes on a fixed 
value, which we shall denote by \/7. Replacing 
(j, k) by (x, y) in (18) we get the parabola 


| n e\* 
(«-a)*="y40)+() , (19) 
16 32 


which is the locus of all points (x, y) for which 
the minimum speed is the same. The curve (19) 
is concave up with negative slope to the left of 
x=a. As we vary 7 we obtain a one-parameter 


) 


a GM PO, ¥) 





9* (ab) 








Fic. 14. 
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family F of parabolas. In particular, if 7=0, we 
have the vertical line x=a through D*. As 7 is 
increased the corresponding parabola is lowered, 
in the sense that the parabola (xSa) for n= 
lies above the parabola for »=n2 when 72> 1. 
(See Fig. 14.) It follows that the point P which 
must be occupied by C, at the starting instant /, 
where the speed is a minimum is the point on 
the circle ! where a member of the family F is 
tangent to T’. We let (x1, y:) designate the 
coordinates of P. The normal to the curve of F 
through P is given by 

y = nx/32(a—x)). (20) 


Setting (x, y)=(x1, y1) in (19) and (20), solving 
for n, and substituting for » in (19), we obtain 


x17(x1—a@) = —2xyyi(yi— 9) +yP(x1—a@). (21) 
The point P is determined by (21) and 
x?+y?=h’. (22) 


The solution of (21) and (22) for (x;, yi) can be 
reduced to that of a sextic in x;. Setting (j, k) in 
(18) equal to (x:, y:) we obtain the minimum 
speed. The minimum kinetic energy which must 
be generated in jumping from A to D* is thus 


E=16m[((b+y1)?+(¢e—x,)*)'— (b+ y;) J. (23) 


The angle a of departure for the most efficient 
jump is given by (17) with (j, k) =(x, y:). If a 
is large enough, e.g., when (a, 6) = (15, 30), the 
energy E of (23) is near the minimum kinetic 
energy Ey which must be furnished by M to 
reach D*. 


3. Combination Fall and Jump 


By falling forward while jumping the body M 
can decrease the kinetic energy which he must 
furnish to reach D*. We shall find a number E,; 
which is numerically greater than this decrease. 

From (3) with g=32 we find that the x-com- 
ponent v, of the velocity v; of M after falling 
forward through an angle @ while standing on 
‘A is 

8Vh 
v,=—— cos 6(1—cos 6)!. 
T 
Setting dv,/d@=0, we find that cos @= 3%, when 
v, is a maximum, whence this maximum com- 
ponent is 
v,= 16(3h)'/9r. 
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Since, in jumping, the angle a which the velocity 
v of departure makes with the x axis is positive, 
the contribution which v; makes to v, is less than 
v,, inasmuch as the component of v; on v is less 
than v, in magnitude. Thus the kinetic energy 
furnished by a fall does not exceed 


¢,=128mh/27 7. (24) 


It follows-from (23) and (24) that the kinetic 
energy which must be furnished by M is at 
least E— Ey. 

In experiments track athletes, aged 20, were 
not able to develop 100 foot pounds of kinetic 
energy after loss of balance. This was determined 
from the maximum range that could be attained 
by athletes jumping after loss of balance. If M 
is a normal individual and the trajectory of M 
is such that E— E,;>100, one may conclude that 
this trajectory was made by a premeditated 
jump. 


4. Steps 


In making a step forward the velocity is 
practically horizontal. Using (j, k) = (0, 0), a=0, 
in (4) we obtain a trajectory with the range 


9) 


c 


x= Vb (25) 


4 


on the line y= —bd. For 3 and 4 m.p.h. (4.4 and 
5.9 ft. per sec.) (25) yields 1.1\/b and 1.5/3, 
respectively. Since the center of gravity C of M 
for a normal body M is at (j, k) =(1.5, 3) in- 
stead of (0,0) at the instant of departure the 
range on y= —d is actually 

lv | , : 

x=1.5+—(b+3)}. (26) 

4 
By comparing x of (26) with the actual range a 
one can determine whether or not the actual 
trajectory to D* could have been attained by a 
step with velocity v. We remark that the first 
step of a walk at v ft. per sec. is much more 
difficult to make than succeeding steps. 


A SPECIAL CASE 


In a New York Life Insurance case, a man VM 
weighing 150 Ib., 5 ft. 8 in. tall, left a ledge A, 
was found on the concrete pavement 30 ft. 
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below A, and with all points of M in the final 
state over 14 ft. from the vertical A Y through A. 
On contact with the pavement the center of 
gravity C of M was traveling nearly along a 
vertical line and the point was at least 14 feet 
from the line AY on contact of M with the 
ground, for otherwise either M would have 
fallen flat on the ground, or would have struck 
on the head or legs, and due to lever action part 
of M would have been found nearer to AY 
than 14 ft. Actually from the crushing state of 
the head it was established that M struck the 
pavement on his head. In analyzing the motion 
we compare ranges on a line y= — 27, supposing 
that the range on this line was 14 ft., so that 
D*=(14, —27). For M, we have h=10/3, r=4. 
By (3) the velocity of departure for the most 
efficient fall had the magnitude 7.9 ft. per sec., 
and by (5) the center of mass for such a fall 


would have been 7.8 feet out on y= —27. This 
means that M would have been found no more 
than 7. ft. out (allowing for end corrections). 
Thus M produced energy of motion himself. 
From (26), with allowances for the ends, M 
would have been found about 7 or 9 ft. out if M 
had stepped forward at 3 or 4 m.p.h., respec- 
tively. Thus M could not have made such a 
step. From (21) and (22), the point (x, y:) 
from which C could have departed with the 
least kinetic energy and reach (14, —27) was 
(3.25, 0.55). The minimum possible energy E 
was over 150 ft. lb., whereas the possible energy 
E; due to falling forward was less than 50 ft. Ib. 
Actually, M furnished at least 150 ft. Ib. himself. 
Because of M’s age and condition, M could 
have furnished only an energy much less than 
100 ft. lb. after loss of balance. Thus M made a 
premeditated jump of considerable magnitude. 





New Instrument Booklets 








Research Progress for June, 1942, symmarizes in an 
article entitled “Inside Stories Told by Atoms” inves- 
tigations in various fields by means of radioactive tracers, 
emphasizing work on metals now being carried out at 
Westinghouse. This is a bi-monthly publication of the 
Westinghouse Research Laboratories, East Pittsburgh, 
Pennsylvania. 


General Radio Experimenter, May, 1942 issue, presents 
the following articles: ‘The Type 727-A Vacuum Tube 
Voltmeter’ and “Dielectric Strength Tests with the 
Variac.” Issued by the General Radio Company, 30 State 
Street,.Cambridge, Massachusetts. 


Cenco News Chats for May, 1942, occasional publication 
of the Central Scientific Company, 1700 Irving Park Road, 
Chicago, Illinois, contains the usual descriptions of ap- 
paratus and a number of interesting photos of the new 
Northwestern Technological Institute at Northwestern 
University. Lord Kelvin is the subject of the cover picture. 
The leading article is entitled ‘‘This is a War of Science.” 


Bakelite Review, April, 1942 issue, presents in a 24-page 
booklet the latest developments in plastics and their 
applications, including Part 2 of a series of three articles 
on ‘‘Designing Molded Plastic Parts.’”’ ‘A Useful Table of 
Plywood Adhesives—Their Characteristics and Methods 
of Handling”’ lists 24 properties of the ten principal types 
of wood adhesives used commercially at present. Quarterly 
publication of the Bakelite Corporation, Unit of Union 
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Carbide and Carbon Corporation, 30 East 42 Street, New 
York, New York. 


Science and Appliance for May, 1942, publication of the 
Ohio State University Research Foundation, Columbus, 
Ohio, features a discussion of the sources of the important 
war metal, magnesium, and its production from dolomite 
and ferrosilicon. 


The 32-page Spring Edition of Inco, issued by the Inter- 
national Nickel Company, 67 Wall Street, New York, 
New York, discusses the salvage of scrap; the part which 
nickel and nickel alloys play in the production of essential 
wartime material, including tungsten and plastics used in 
the construction of military and naval airplanes; and 
describes the development of a new type of rivet which 
eliminates blind fastening. 


The Spring, 1942, edition of Modern Precision has just 
been issued by Leeds & Northrup Company, 4934 Stenton 
Avenue, Philadelphia, Pennsylvania. This is a 16-page paper 
containing brief news items about L & N instruments, 
controls, and heat-treating methods. It supplements the 
regular catalogs, bulletins, and circulars which will con- 
tinue to be issued as before. 


In the Tin Research Institute’s Report for 1941 it is 
stated that since Japan began the Far Eastern campaign 
the main problem has become one of how to make the best 
use of the supplies available by curtailing the use of tin in 
commodities not essential to the war effort; by improving 
the processes in essential uses; and by recovering tin. The 
Institute has been appointed by the Ministry of Supply 
to advise on all technical questions affecting the use of tin. 
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Collaboration with manufacturers and with government 
departments in the application of electro-tinning for parts 
of armaments was continued during the year, and there has 
also been considerable interest in electro-tinning for small 
soldering tags, terminals, and clips. Other sections of the 
Report review the progress made in research on tin-rich 
bearing alloys, foil, and bronzes. Copies may be obtained 
free of charge from the Tin Research Institute, Fraser 
Road, Greenford, Middlesex, England. 


Low-voltage air circuit breakers designed for use on 
alternating current circuits up to 600 volts and direct 
current circuits of 250 or 750 volts are described and. illus- 
trated in a new 12-page publication, Catalog 2150, just 
released by the Roller-Smith Company, Bethelehem, Penn- 
sylvania. The catalog lists closing and tripping currents, 
shows typical methods of making connections to large 
breakers, gives dimensions of both manual and electrical 
designs, and indicates the number and size of terminals for 
various capacity breakers. 





Innovations in Instruments 








Microfilm Reader 


The Spencer Microfilm Reader, developed as a result of 
studies by the Committee on Scientific Aids to Learning, 
is the first inexpensive instrument designed to reproduce 
the almost microscopic characters of microfilm with bril- 
liant fidelity. The instrument consists of a projection head, 
a glass film book, and a shadow box and screen. The pro- 
jection head incorporates an inexpensive spotlight type of 
longlife bulb, the proper condensing system, and a specially 
designed projection lens. This head is held by spring clips 
to the shadow box, which also serves as a support, and may 
be rotated to project the film in either a vertical or hori- 
zontal meridian as desired. Operation of the instrument is 
simple. The microfilm is placed in a convenient glass film 











book which consists of two pieces of glass, hinged and 
bound at the edges with transparent Cellophane tape. This 
tape serves to space the film properly to protect it from 
scratching and to orient the film in relation to the focal 
plane of the objective. The image is projected onto the 
screen and is so protected by the shadow box that a 
darkened room is unnecessary. The image is flat and undis- 
torted and should not cause optical fatigue or discomfort. 

The instrument is regularly supplied with one bulb and 
one glass film book, at a list price of $37.50. A roll film 
attachment designed to accommodate rolls of film up to 
100 feet in length is listed at $15.00. Further details may 
be obtained from the Spencer Lens Company, Buffalo, 
New York. 


New Cellulose Product 


A new cellulose product, Onco “V’’, has been developed 
by Brown Company, 500 Fifth Avenue, New York, New 
York. The cellulose base for Onco ‘‘V”’ was originally dis- 
covered during research experiments on filter materials for 
gas masks. First commercial use of Onco ‘‘V” was as an 
insole fabric to replace latex-impregnated materials com- 
monly used before the rubber restrictions. It is announced 
that this product can be used as a substitute for leather 
and cork in the manufacture of luggage and linoleum, 
gaskets, and leather goods. 


Improved Blackout Bulb Unit 


Changes in blackout bulb specifications are announced 
by the Wabash Appliance Corporation, Brooklyn, New 
York, whose silver-lined blackout bulb, placed on the 
market in early January, has been tested in actual city- 
wide blackouts in practically every. state. 

The most important specification change is in color of 
light from blue and red to the deep orange recommended 
by the Office of Civilian Defense. The new bulb is smaller 
in size, with a reduced current consumption of 15 watts. 
The built-in reflector has been eliminated, and an improved 
type of heavy black silicate coating is used to prevent light 
leakage. The bulb will fit any household socket and will 
list at 45¢. 


New Wood Table Radios 


Substantial savings of formaldehyde, resins, and other 
chemicals and materials used to make plastics needed for 
national defense have been made possible in the production 
of two new RCA Victor wood table radios to replace several 
plastic-housed instruments of similar type. The instru- 
ments, designated as Models 35X and 36X, operate on 
a.c.-d.c. current, and both have Underwriters’ approval. 
Standard broadcasts and police calls are received. 

Model 35X is housed in a solid walnut cabinet and has 
a new Spread-Vision dial which, with extra large tuning 
knobs and 12:1 ratio vernier tuning, provides precision 
tuning. The 36X is similar to the famed RCA Victor Super 
Six instrument in a wood cabinet. It also employs the new 
Spread-Vision dial. A 6-tube instrument, the 36X has a 
2-point high frequency tone control and an extra large 
built-in Magic Loop antenna. Both new models employ 
Magnetite core transformers for greater selectivity. 
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